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1.  INTRODUCTION 


This  final  technical  report  is  concerned  with  work  done  under  Grant 
AFOSR  78-3635,  A  Study  of  Hybrid  Computing  Techniques  for  Transonic  Flow 
Fields.  The  grant  was  started  on  August  1,  1978  and  was  to  run  for  a  period 
of  14  months.  A  three  month  no-cost  extension  was  granted  so  that  the  final 
date  for  technical  work  was  December  31,  1980. 

The  purpose  of  this  work  is  the  investigation  of  the  use  of  asymptotic 
techniques  to  derive  analytical  solutions  for  use  in  conjunction  with  numerical 
methods,  to  improve  the  accuracy  of  solutions  or  to  reduce  computing  times  or 
both.  The  two  problems  chosen  for  study  are 

(1)  The  use  of  asymptotic  far  field  solutions  in  viscous  flow  field 
computations 

(2)  The  accurate  determination  of  the  location  of  a  shock  wave  in  a  flow 
field,  using  composite  asymptotic  solutions  in  the  neighborhood  of 
the  shock  wave. 

The  first  problem  is  concerned  with  the  use  of  analytical  far  field  solutions 
rather  than  boundary  conditions  at  infinity  to  set  conditions  which  the  numerical 
computations  must  approach  as  the  solution  proceeds  away  from  the  body  in 
question.  For  inviscid  flow  fields,  such  outer  solutions  have  been  used  in  the 
past  with  considerable  success;  as  points  farther  from  the  body  are  considered, 
the  numerical  solution  is  compared  with  the  far  field  solution  and  computation 
may  be  stopped  when  the  difference  between  the  two  is  within  the  allowable  error. 

In  the  present  work,  covered  in  Section  2,  far  field  solutions  for  viscous  flow 
fields  are  considered  for  axis  ymme  trie,  super  sonic  and  transonic  flows.  For  a 
given  set  of  parameters,  the  results  allow  determination  of  the  length  scales  at 
which  viscous  effects  are  important  and  therefore  show  the  general  nature  of  the 
solutions  valid  at  the  various  length  scales.  Then,  depending  upon  the  accuracy 
desired,  a  decision  can  be  made  as  to  the  type  of  outer  solution  to  be  used. 
Analytical  solutions  are  obtained  for  the  case  judged  to  be  probably  the  most 
important. 


In  the  second  problem,  covered  in  Section  3,  the  flow  problem  chosen 
as  an  example  problem  is  that  of  steady  transonic  flow  through  a  two-dimensional 
channel  in  which  a  shock  wave  is  located.  Because  several  regions  of  flow  must 
be  considered,  adding  to  the  general  complexity  of  the  solution,  the  solutions  in 
this  section  are  presented  in  detail. 

Because  two  independent  problems  were  studied,  the  sections  in  which 
each  is  considered,  Sections  2  and  3,  are  independent  and  self-contained. 

The  work  described  herein  represents  only  the  initial  phases  of  a  study 
of  hybrid  computing  techniques  in  that  only  analytical  solutions  are  considered 
for  only  two  problems. 


2.  FAR -FIELD  DIFFUSION  EFFECTS 
IN  AXISYMMETRIC  SUPERSONIC  AND  TRANSONIC  FLOW 

A.  F.  Messiter 

Introduction 

The  propagation  of  a  weak  pressure  disturbance  in  a  gas  which  initially 
is  at  rest  and  has  uniform  properties  is  described  by  the  linear  wave  equa¬ 
tion,  provided  that  the  distance  through  which  the  wave  travels  is  not  too 
large*  Nonlinear  and  diffusion  effects,  however,  may  eventually  become  im¬ 
portant;  the  details  differ  for  plane  and  cylindrically  symmetric  flows.  The 
present  work  is  concerned  with  axis ymme trie  supersonic  and  transonic  flow 
past  a  slender  body  of  revolution  in  the  case  when  nonlinearity  is  an  essen¬ 
tial  feature  at  distances  smaller  than  the  distance  at  which  a  fully  viscous 
region  appears. 

The  inviscid  far  field,  nonlinear  because  shock  waves  and  characteris¬ 
tics  are  slightly  inclined  from  their  linearized  positions,  was  studied  by 
Whitham^;  a  summary  appears  in  his  book.  ^  A  detailed  account  of  vis¬ 
cosity  effects  on  sound  waves  was  given  by  Lighthill.  ^  In  particular,  he, 

(4) 

and  later  Ryzhov  and  others,  noted  that  for  cylindrically  symmetric 
motions  a  fully  viscous  region  finally  develops,  which  is  described  by  an  equa¬ 
tion  similar  to  the  Burgers  equation  but  modified  because  of  the  geometry. 

The  velocity  field  beyond  this  fully  viscous  region  was  studied  by 

^Whitham,  G.  B.,  The  Flow  Pattern  of  a  Supersonic  Projectile,  Comm. 

Pure  Appl.  Math.,  5_  (1952),  301-348. 

^Whitham,  G.  B. ,  Linear  and  Nonlinear  Waves.  John  Wiley  and  Sons,  1974. 
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Lighthill,  M.  J.,  Viscosity  Effects  in  Sound  Waves  of  Finite  Amplitude, 
in  Surveys  in  Mechanics,  eds.  G.  K.  Batchelor  and  R.  M.  Davies,  Cam¬ 
bridge  University  Press,  1956,  pp.  250-351. 

a 

Ryzhov,  O.  S.,  Influence  of  Viscosity  and  Thermal  Conductivity  on  Propa¬ 
gation  of  Sound  Impulses,  P.  M.  M. ,  30^  (1966),  296-302.  Engl,  transl. 

J.  Appl.  Math.  Mech. ,  32.(1966),  362-369. 


S  ala  the  ,  Chong  and  Sir  ovich,  and  Sanchez -Palencia -Hubert.  Chong 

(8) 

and  Sirovich1  7  also  considered  the  combined  effects  of  diffusion  and  non¬ 
linearity.  For  a  body  moving  at  sonic  speed,  these  effects  have  been  studied 

by  Szaniawski^  and  Ryzhov  and  Shefter(^  Review  articles  by  Sichel^^ 

(12) 

and  by  Ryzhov  include  discussions  of  these  flows  and  list  additional 
references. 

The  present  study  has  two  main  purposes.  First,  for  small  but  nonzero 
viscosity  and  thermal  conductivity,  the  far -field  solution  for  inviscid  super¬ 
sonic  flow  can  be  modified  to  include  the  details  of  the  velocity  variation 
through  the  front  and  rear  shock  waves.  This  requires  calculation  of  second  - 
order  solutions  for  the  velocity  in  the  region  between  the  shock  waves  and  for 
the  velocity  distribution  within  the  shock  waves,  so  that  the  uncertainty  in 

shock -wave  location  is  smaller  than  the  shock-wave  thickness.  These 


Salathe,  E.  P.,  The  Fundamental  Matrix  in  Three-Dimensional  Dissipative 
Gasdynamics,  J.  Fluid  Mech. ,  39^(1969),  209-225. 

6 

Chong,  T.  H.  and  Sirovich,  L.  ,  Structure  of  Three-Dimensional  Supersonic 
and  Hypersonic  Flow,  Phys.  Fluids,  JL_3^  (1970),  1990-1999. 

7 

Sanchez -Palencia -Hubert,  J. ,  Asymptotic  Study  at  Infinity  of  Supersonic 
Flows  of  a  Dissipative  Fluid,  Int.  J.  Eng.  Sci. ,  1^(1976),  567-584. 

8 

Chong,  T.  H.  and  Sirovich,  L.. ,  Nonlinear  Effects  in  Steady  Supersonic 
Dissipative  Gasdynamics.  Part  2.  Three-Dimensional  Axisymmetric  Flow, 

J.  Fluid  Mech.,  58.(1973),  53-63. 

9 

Szaniawski,  A.,  The  Asymptotic  Structure  of  Weak  Shock  Waves  in  Flows 
over  Symmetrical  Bodies  at  Mach  Number  Unity,  Acta  Mech.,  (1968), 
189-203. 

^Ryzhov,  O.  S.  and  Shefter,  G.  M. ,  On  the  Effect  of  Viscosity  and  Thermal 
Conductivity  on  the  Structure  of  Compressible  Flows,  P.M.M.,  2jB  (1964), 
996-1007.  Engl,  transl.  J.  Appl.  Math.  Mech.,  28  (1964),  1206-1218. 

^Sicbel,  M. ,  Two-Dimensional  Shock  Structure  in  Transonic  and  Hypersonic 
Flow,  in  Advances  in  Applied  Mechanics.  11,  ed.  C.  S.  Yih,  Academic  Press, 
1971,  pp.  131-207. 

^Ryzhov,  O.  S.,  Viscous  Transonic  Flows,  in  Ann.  Rev.  Fluid  Mech. ,  10 
(1978),  65-92. 
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solutions  are  derived,  the  second -order  shock -wave  position  is  obtained, 
and  a  simple  multiplicative  composition  then  gives  a  solution  which  is  uni¬ 
formly  valid  in  the  streamwise  coordinate.  Second,  the  flow  properties 
depend  on  the  Mach  number,  the  Reynolds  number,  and  the  body  slender¬ 
ness  ratio,  and  different  approximations  are  needed  for  different  parame¬ 
ter  and  coordinate  ranges.  In  particular,  at  distances  where  nonlinear 
effects  first  become  essential  the  shock  waves  are  thin  and  one -dimension¬ 
al;  somewhat  closer  to  the  body  the  shock -wave  structure  becomes  two- 
dimensional  and  somewhat  farther  away  the  shock  waves  can  no  longer 
be  considered  thin.  In  the  existing  literature,  for  both  supersonic  and 
transonic  flows,  the  identification  of  different  approximations  and  their 
domains  of  validity  is  incomplete.  Several  extensions  and  corrections 
are  given  here,  so  that  the  dominant  physical  effects  can  be  identified 
as  functions  of  position  and  of  the  nondimens ional  parameters. 


Basic  Equations 


Dimensional  quantities  are  indicated  by  an  overbar;  the  subscript  oo 

refers  to  values  in  the  undisturbed  uniform  flow  ahead  of  the  body.  Thus, 

e.  g. ,  u  ,  a  ,  p  ,  p.  ,  and  X  are,  respectively,  the  undisturbed  values 
00  oo  oo  OO  00 

of  the  flow  speed,  sound  speed,  fluid  density,  and  first  and  second  viscosity 
coefficients.  The  nondimens  ional  flow  speed  q,  sound  speed  a,  enthalpy  h, 
density  p,  pressure  p,  and  stress  tensor  r  are  defined  by 


P 


P 


P-2 

u 

00  CO 


T 


P  u 
00  oo 


(2.1) 


The  Reynolds  number  Re  is  based  on  the  body  length  L  and  a  reference  vis¬ 
cosity  coefficient  specified  later  to  be  2ix  +  X  ;  the  same  reference  vis- 

00  oo 

cosity  is  used  for  the  Prandtl  number  Pr  and  the  nondimensional  first  and 
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(2 


second  viscosity  coefficients  p.  and  X  .  In  terms  of  these  nondimens ional 
variables,  the  differential  equations  describing  the  fluid  motion  are 

div  p  q  =  0 

pq  •  Vq  +  grad  p  =  Re  ^div  t  (2 

pq*VhQ  =  Re  ^div(Pr  ^grad  h)  +  Re  ^div(Tq)  (2 

where  hQ  is  the  nondimens  ional  total  enthalpy.  Also 

t  =  p.  (def  q  )  +  X  (div  q  )I  (2 

where  def  q  and  I  are  the  deformation  and  identity  tensors.  A  useful  com¬ 
bination  of  the  differential  equations  is 

a  div  q  -  q  •  V  ?"•  =  (y-l)q-Vh  -Re  ^  q  *  div  t  (2 

t  op 

A  perfect  gas  with  constant  specific  heats  is  assumed;  y  is  the  ratio  of 
specific  heats. 

Cylindrical  coordinates  x  =  x/  L,  r  a  r/  L,  and  8  are  introduced,  with 
the  origin  at  the  nose  of  the  body  and  x  measured  in  the  direction  of  the  un¬ 
disturbed  stream.  For  axisymmetric  flow  the  velocity  vector  is 


q  =  (l+u)ex  +  ve  r 

2  2 

where  e  and  e  are  unit  vectors;  if  the  body  is  slender,  u  +  v  <<  1. 
x  r 


(2.7) 


Again  for  axisymmetric  flow 


div  t  =  ( 


1  8,rTxr>  - 


+  r  at 


ax  ,  a(rT  ) 

,  ,  rx  I_  rr 

ex  3x  r  3r 


99  - 
~)er 


(2.8) 


where,  with  ©  =  div  q  , 


t  =  “  +  \© 

xx  9x 


t  =  Z\i  ~  +  \© 
rr  3r 


.  3v  3u 

Txr  =  Trx  -  111  to  *  57  > 


t99  *  2*7  + 


(2.9) 
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Supersonic  Far  Field  with  Small  Diffusion 


For  r  =  0(1)  the  flow  past  a  slender  body  is  described  in  a  first  approx 
imation  by  a  velocity  potential  which  satisfies  the  wave  equation  in  the  coor¬ 
dinates  x  and  r.  For  a  body  of  revolution  at  zero  incidence  the  flow  is  axi- 

symmetric,  and  the  solution  is  expressed  in  terms  of  a  distribution  of  sour- 

2 

ces  along  the  axis,  with  strength  equal  to  the  rate  of  change  €  S'  (x)  of  the 
body  cross-section  area,  where  €  «  1.  As  r  oo,  this  solution  gives 

2 

~  J0 


u  ~  - 


2ir(2  Br) 


J 


x  -B  r  _  *  * 

SM(x  )dx 


(x-Br  -x 


*,1/2 


(2. 


for  x  >  Br,  where  B  =  (M^  -  1)^  ^  and  M  =  u  /  a  ;  the  body  has  been 

oo  oo  oo  oo 

taken  to  have  a  pointed  nose  so  that  S1  (0)  =  0.  The  result,  however,  is  not 
uniformly  valid  for  large  r  because  nonlinear  and  viscous  effects  have  not 
been  taken  into  account. 

In  the  case  considered  here,  the  first  error  encountered  as  r  becomes 
large  arises  because  the  characteristics  have  slopes  slightly  different  from 
the  undisturbed  value  dx/ dr  =  B.  The  weak  shock  waves  from  the  nose  and 
rear  of  the  body  likewise  have  slopes  slightly  different  from  B.  At  a  still 
larger  distance,  viscous  forces  would  also  have  to  be  included  in  the  first 
approximation.  If  the  displacement  of  the  characteristics  from  their  linear¬ 
ized  locations  is  0(1)  when  1/ r  =  0(5),  where  6  =  6(c),  it  is  convenient  to 
introduce  a  coordinate  r}  =  6(  e)r.  Since  x  -  Br  =  0(1)  in  the  disturbed 
portion  of  the  flow,  a  second  suitable  coordinate  is  £  =  x  -  Br.  The  equa¬ 
tions  are  then  studied  in  a  limit  as  €  0  with  the  coordinates  £  =  x  -  Br 

and  rj  =  5(e)r  held  fixed.  As  r\  0,  the  solution  found  in  terms  of  |  and  rj 

must  match  with  the  solution  (2. 10)  which  is  correct  for  1  «  r  «  6 

2  -l/  2 

Since  the  solution  (2.10)  gives  u  =  0(e  r  ),  it  follows  from  the  matching 
that  u  =  0(e^6^^)  when  r  =  0(6  *).  The  requirement  that  a  nonlinear 

4 

term  appear  in  the  differential  equation  for  the  potential  then  gives  5=  0(e  ); 

4 

the  simplest  choice  is  5(e)  =  e  .  It  is  also  seen  from  this  equation  that 

-1  8 

the  next  terms  in  u  are  0(Re  )  and  0(e).  In  the  case  considered  here. 


10) 
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6^  «  Re  ^  that  is,  a  limit  is  taken  with  €  -►  0  and  Re  oo  such  that 

4  8 

€  Re  oo  and  e  Re  -*■  0.  The  coordinates  and  the  expansion  for  the  velocity 


are  then  given  by 


I  =  x-Br, 


T)  =  €  r 


(2. 11) 


4  -1 

u  =  e  ,  n)  +  Re  u2(£  ,17)  +  — 

The  shock  wave  from  the  nose  of  the  body  is  located  at  %  =  £  g(r?)»  where 

§g  =  !sl(rj)  +  (€4Re)'1eg2(n)  +  ... 


(2.12) 


(2.13) 


The  rear  shock  wave  is  represented  in  a  similar  way.  Some  general  features 

4 

of  the  flow  for  e  r  =  0(1),  and  in  the  more  distant  fully  viscous  region  dis¬ 
cussed  below,  are  shown  in  Fig.  1;  coordinates  used  here  and  later  are 
shown  in  Fig.  Z. 

The  differential  equation  for  u^,  found  from  Eq.  (Z.  6)  in  the  manner 
described  above,  is 


U  (v+i)*T 

2%  +  “  +  ~  Vii 


(2. 14) 


1/  Z 

After  multiplication  by  77,  the  differential  equation  states  that  77  is 

constant  along  characteristics  d£/dr?^^  =  (y+1)M^B  ^u^.  The  solution 

which  satisfies  the  matching  condition  for  77  0  is  given  by 


Uj  =  -(2Brj)'l/2F(X), 


X  =  4  +kr?l//2F(X) 


™  -  h  J* 

0  (X-x)' 


(2.15,16) 


(2.17) 


where  k  =  (y+l)M4  /(2B2)^2.  This  is  exactly  Whitham' s  solution.  2 ^  The 

GO 

shape  X  =  X  (17)  of  the  front  shock  wave  is  most  easily  found  by  use  of  the 
s 

known  results  that  a  weak  shock  wave  approximately  bisects  the  angle  be¬ 
tween  upstream  and  downstream  characteristics;  from  Eqs.  (Z.  13)  and  (Z.  16), 


it  follows  that 


*;!<•»  *  -  l-k,-1/2F(Xs) 


(2.18) 
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The  result  found  for  X  (77 )  is  given  implicitly  by 

s 

x 

7ki)1//2F2(X)  =  j  F(X*)dX*  (2.19) 

2  s  0 


Evaluation  of  u.  near  the  shock  wave  requires  an  expansion  in  the  form 


s2 


Ujd,!))  =  Ul^9l'T7)+~4 -  Ul|  ^  sT  ^  +  ‘  '  *  +  ^  3^4  ^  si'  ^  + 

€  Re 


(2*  20) 


where  u^  (£  77)  is  found  by  differentiating  the  solution  (2.15)  for  u^  and 

the  definition  (2.16)  of  X,  and  combining  with  the  derivative  of  Eq.  (2.18) 
for  the  shock-wave  shape;  the  result  is 


4B  d_ 

<V«)M4 

00 


In  F(X  ) 
s 


(2.21) 


For  the  calculation  of  u^,  the  first  approximation  to  the  diffusion 

terms  is  obtained  using  - - Bu  and  3u/3r - B3u/3x.  From  Eq.  (2.  4) 

4  -1 

it  is  found  that  changes  in  h^  along  a  streamline  are  0(e  Re  );  entropy 
changes  are  also  of  this  order.  It  follows  from  the  vorticity  equation, 

g 

found  by  taking  the  curl  of  Eq.  (2.  3),  that  the  vorticity  is  0(e  /Re),  and 
so  the  flow  remains  irrotational  to  the  order  required  here.  The  differen¬ 
tial  equation  for  u^  is 


(y+l)M^ 

00  , 

~  (V2>5 


Pr 


)MLU4! 


(2.22) 


In  the  right-hand  side,  the  viscosity  coefficients  appear  in  the  combination 

Z\i  +  X  ~  1,  the  reference  value  having  been  taken  to  be  Z\x  +  X  .  Equa¬ 
l/2  °°  °° 

tion  (2.22)  can  be  rewritten  as  an  equation  for  ry  u^  in  terms  of  coor¬ 
dinates  ry*^  and  X,  where  X  is  constant  along  a  characteristic.  The  solu¬ 
tion  can  then  be  found  as 


u 


2 


-  (1  + 


ni) 

Pr  ' 


M4  F" 

_ 00 _ 

B(2Bn)1//2(l-kn1//2F'  )k2F’  2 


ln(l 


■kn1/2F' 


+ 


*  n1/2F'  ) 

i-kn1//2F’J 


(2.23) 


where  F'  =  F1  (X);  an  arbitrary  function  of  X  has  been  chosen  as  that 
remains  bounded  as  r)  0* 


Shock* Wave  Structure  and  Location 


In  order  to  find  the  second  term  £^(77)  in  the  shock-wave  shape,  it  is 

necessary  to  obtain  both  first-  and  second -order  solutions  for  the  velocity 

distribution  within  the  shock  wave.  Coordinates  x  and  y  are  defined  to  be 

-4 

of  order  one  at  points  within  the  shock  wave  when  r  =  0(e  )  and  are  meas¬ 

ured,  respectively,  normal  to  and  along  a  curve  which  remains  within  the 
shock  wave  as  e  0  and  Re  —  oo  (Fig.  2).  This  curve  is  chosen  such  that 
the  normal  velocity  component  is  sonic  at  every  point,  and  its  angle  of  in¬ 
clination  g(rj)  from  the  radial  direction  is  related  to  the  expansion  (2.15)  of 

4  s(n)  by 


ctn  p(rj)  =  B  +  e4|  ^(n)  +  Re  1|  ^2(rj)  +  •  •  •  (2.24) 

The  coordinate  transformation  from  x,  r  to  x,  y  can  be  carried  out  in  three 

steps.  First,  x  and  r  are  replaced  by  £  -£  ^ ( 77 )  and  17;  next,  at  each  point 

the  coordinate  directions  are  rotated  through  the  angle  g(n);  and  finally, 

the  coordinate  measured  normal  to  the  shock  wave  is  stretched  by  a  factor 

A  (y)  «  1,  where  A  (y)  is  a  suitably  defined  shock-wave  thickness,  as  ob- 
s  s 

tained  below.  The  transformation  of  derivatives  can  be  summarized  in 
the  form 


8_ 

3x 


~  dA 

1.8.4  ,  a  x  s  a  , 

~  sin  r  TZ  +  €  cos  b  (  —  -  7 —  ~ —  rr: ) 
A  p  9x  p  8y  A  dy  9x 

s  s 


(2.25) 


a_ 

ar 


1  a  .  4  .  .  a  x  ^s  a  . 

-  T*  cos  b  7—  +  €  sin  g  (  — -  -  - — - - — ) 

A  p  8x  p  3y  A  dy  85? 

S  9 


(2.26) 


The  velocity  perturbations  U  and  V,  in  the  x  and  y  directions  respectively, 
are 


10 


f 


U  =  u  sin  g  -  v  cos  g  i 

V  =  u  cos  g  +  v  si  n  g 

g 

Changes  in  entropy  and  in  total  enthalpy  within  the  shock  wave  are  0(e  ). 

The  vorticity  is  extremely  small,  of  order  since  the  shock-wave 

-1  4 

curvature  is  very  small,  and  it  follows  that  A  V-*  ~  e  U  .  The  differen- 
7  s  5?  y 

tial  equation  for  U  is  found  to  be 


ry  ^  •  b2)  r  - w+1)  t  uus 

M  s  s 

00 


(2.27) 

(2.28) 


4_  |  TT  ~  dA 

J-+2MU  -M  “  — 3  uJ 
M2  |  T?  00  y  oo  Ag  dy  x 

oo  L 


(2.29) 


where  dr^/dy  =  M^.  Equation  (2.29)  is  to  be  studied  in  the  limit  as  e 
and  Re  oo  with  coordinates  2  and  y  held  fixed. 

The  velocity  U  is  represented  in  the  form 


-  1  +  CJ  (52)  -  —  +  — 

vi  e  R«6'  (r?)  €  Re 

o  X 


U2(x;y)  +  .  .  . 


(2.30) 


where  2AU  is  the  first  approximation  to  the  jump  in  U  across  the  shock  wave, 

and  |  1  ( T7 )  is  given  by  Eq.  (2.18).  The  shock-wave  thickness  A  used  in  the 
sl  s 

definition  of  x  is  found  by  equating  the  coefficients  of  the  nonlinear  and  diffu¬ 
sion  terms  in  Eq.  (2.29),  with  a  factor  of  two  included  for  later  convenience. 

Thus  A  and  A  are  given  by 
u  s  ° 


A  M  FIX) 

4  oo  s 

2{2Bn)1/2 


2 

(■y+l)A 


Pr  f  Re 


(2.31) 


As  for  a  one^dimensional  shock  wave  of  the  same  strength,  e  A  is  the  first 

u 

approximation  to  the  difference  between  the  normal  velocity  component 
ahead  of  the  shock  wave  and  the  critical  sound  speed  based  on  that  velocity. 


11 


Since  x  has  been  defined  such  that  the  normal  velocity  component  is  sonic  at 
5?  =  0,  it  follows  that  U^(0)  =  0;  similarly,  for  the  second  approximation  it 
is  seen  that  U^tO)  =  0.  As  x  -  oo,  —  1  and  | 


The  differential  equations  for  and  are 


-  2U  U  =  0 
125?  1  12 


(2.  32) 


dA  .  dA 

U2SR  - 2«W*  * A ; f  «■*  +  vi  *  sf  “i* >  (2-331 


-1  -2  8  -1 

where  A  =  2B(v+l)  M  €  ReA  A  =  0(1).  The  solutions  which  satisfy 
*  OD  S  U 

the  required  boundary  conditions  are  given  by 


UL(x)  =  -  tanh  x 


(2.34) 


“(4  cosh22)U2(x;n) 


(e2*  -  2x  -  1)~  In  F(X  ) 
arj  s 


-  (e^X  +  x  -  l)x  ^  ln[t)  ^^FfX^)] 


-  {  (sinh  22)ln(l  +  e2^)  +  J  t  Hn(l +t)dt}—ln[n^  F(X^)j 


A5^(n) 


(x  +  sinh  x  cosh  x) 


(2.35) 


The  result  for  U  (5c),  with  the  definition  of  A  given  above,  is  the  familiar 

1  3 

solution  for  a  weak  one -dimensional  shock  wave.  As  x  -*•  oo,  the  solution 
for  U^fx)  has  the  form 

d  ^s2(r,) 

a.  36) 

sx 

After  the  coordinate  transformation  is  taken  into  account,  the  first  term  in 
U,  is  found  to  match  correctly  with  the  result  of  Eq.  (2.21).  Matching  the 
constant  terms  of  order  Re  in  the  expansions  (2. 12)  and  (2.  30)  of  u, 

evaluated  for  §  —  §  and  for  x  —  oo  respectively,  involves  terms  which  appear 

s 


in  Eqs.  (2.20),  (2.23),  (2.  30),  and  (2.36).  The  condition  obtained  gives  a 
first-order  differential  equation  for  the  second  term  £  ?(rj)  appearing  in 
the  shock-wave  shape  (2.13): 


M4  F 
00 


(1+ pT) - 2L77I — 2  If  +fL(1  -k^1/2F’) 

Pr  2B(l-kr71/2F')2  2  F 


(1  -kn^F'  )F" 
2kn1/2  F'2 


ln(l  -  k r)1//2F<  ) 


(2.37) 


where  F,  F’  ,  and  FM  are  evaluated  at  the  shock  wave  X  =  X  (17).  As  t\  0, 

1/2  s 

F(X  )  =  0(X  )  and  X  =  O(rj),  and  the  right-hand  side  of  Eq.  (2.37)  is 

-1/2  s  3 

0(rj  ).  Integration  of  Eq.  (2.37)  should  be  carried  out  from  n  =  0,  so 
that  £  ^(rj)  13  bonded  as  n  0. 

A  solution  for  the  velocity  u  which  is  uniformly  valid  for  all  %  can  be 

4 

constructed  as  the  product  of  e  with  suitable  factors  which  account  for 
the  front  and  rear  shock  waves.  For  the  rear  shock  wave,  a  solution  can  be 
derived  in  a  form  analogous  to  Eq.  (2.  30),  in  terms  of  a  coordinate  Z 
analogous  to  x.  The  composite  solution  then  can  be  written  as 


u  ~  —  c  (1  +  tanh  x)(l  -  tanh  x^)u^(£  ,  77) 


(Z.38) 


To  the  degree  of  approximation  needed  here,  the  coordinate  transformation 
given  by  Eqs.  (2.25)  and  (2.26)  reduces  to 


M  A  (y)x  =  %  -  i  An)  -  — — 

00  3  J  *  Si  '  4_ 

€  Re 


Ss2<n>. 


y  =  M  n 
7  00 


(2.39) 


where,  as  shown  in  Fig.  2,  x  and  y  are  measured,  respectively,  normal  to 

4  -1 

and  along  the  curve  £  =  %  (17)  +  (c  Re)  £  ?(nh  a.  corresponding  definition 

^  si 

can  be  given  for  x^.  Since  |  ^  has  been  calculated,  the  error  in  the  location 

of  the  curve  §  =  §  (r])>  and  therefore  in  the  origin  for  x,  is  of  higher  order 

3  -4  -1 

than  the  shock-wave  thickness  A  =  0(c  Re  ). 

s 


Diffusion  Effects  Elsewhere 

Small  values  of  X  identify-  characteristics  which  intersect  the  body  near 

the  nose,  which  has  been  assumed  pointed  so  that  Sn(x)  remains  finite  as 

x  0.  Thus  for  X  «  1  the  flow  differs  only  slightly  from  the  flow  past  a 

slender  circular  cone  having  surface  r  =  e  X  x.  The  definition  (2.17)  gives 

F(X)  ~2\  V/2  as  X  —0,  and  the  first  approximation  (2.18)  to  the  shock - 

9  2  4 

wave  shape  shows  that  X  (17)  ~jk  X  17  as  17  -*  0.  The  expansion  (2.13)  for 

S  4 

the  shock-wave  shape  becomes 


ftf(n)  =  -  J k2X4r?  -  A(1  +  ^-)  — — —  (-1  +|-ln  3)  +  ..  .  (2.40) 

3(y+1) 

4  -1 

where  A  =  (e  Re)  .  The  second  term  remains  small  in  comparison  with 
the  first  term  provided  that  rj  »  A.  As  rj  0  the  largest  terms  in  the 
differential  equation  (2.29)  become 


tJL\  ?n+JL wJLx  --2B  € 

A  xx  A  A  x  ,  ,....2 

u  u  u  (y+l)M 


A 
co  u 


U/A 

n/A, 


+  2  M 


oo 


3(U/Au) 

a(y/As) 


(2.41) 


Since  both  A  and  A^  approach  constant  values  as  17  -*  0,  the  right-hand  side 

is  small,  and  the  first  approximation  leads  to  the  one-dimensional  velocity- 

distribution  U  ~A  (-1  +  U.),  only  if  n  >>  A.  If  rj  =  0(A),  and  sc 
-8  -L  a 

r  =  0(e  Re  J,  all  the  terms  shown  in  Eq.  (2.  41)  must  be  retained  in  the 

first  approximation.  The  differential  equation  obtained  is  related  to  Burgers' 

equation,  but  contains  an  additional  term  because  the  flow  is  axisymmetric. 

As  n  -*  oo,  X  (rj)  -**  X  =  constant,  where  X  =  X  defines  the  charac- 
so  o 

teristic  at  which  =  0.  Then  Eq.  (2.19)  gives 


X^ 

F2(X  ) - \j-  J  F(X*)dX*  (2.42) 

kr7  '  0 

so  that  F(X  )  =  O(r)  *^).  Also  |  .  =  0(17^),  A  =  O («%  and 

3^4  1  ^ 

A  =  0(An  ),  where  again  A  =  (e^Re)  ,  proportional  to  the  ratio  of  second- 

s 

order  to  first-order  terms  in  the  preceding  results.  It  follows  (Fig.  1)  that 
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the  shock-wave  thickness  is  no  longer  small  in  comparison  with  the  distance 

-2 

between  the  front  and  rear  shock  waves  when  rj  =  0(A  ).  Then  also 

|  s  0(A  and  u  =  0(e  A^2)  in  this  fully  viscous  region.  The  largest 

terms  in  the  differential  equation  (2,6),  all  of  the  same  small  order  of 
magnitude,  give  an  equation 


,  (y+l)M 

4_  .  u  v  .  1  oo 

€  (2u„  +  «  }  + - 5 - 

rj  rj  a 


1  f1  ,  V  “1  V _ CO 

Te(1+  P 7]~  5  { 


(2.43) 


which,  like  Eq.  (2.41),  is  similar  to  Burgers’  equation. 

At  still  larger  distances,  for  rj  »  A  ^ ,  the  linear  terms  are  all  of  the 

same  order  of  magnitude  if  |  =  0(A and  the  nonlinear  terms  are 

-4  -3/2 

smaller  provided  that  €  A  u  «  1.  The  result  is  the  one -dimensional 
diffusion  equation  for  (A^  ?])^^(e^A^^)  \i.  Self-similar  solutions  have 


the  form 


,  4  3/2,  -1  /A  2 
(e  A  )  u  =  (A  i 


where  t  =  A  /rj^^  and  g  satisfies  the  ordinary  differential  equation 


(2.44) 


2  (1  +p7')  b"29  g"  +  F  tg'  +  o;g  =  0 


(2.45) 


The  nonlinear  terms  are  of  higher  order  if  a  >  -  1/2.  If  at  large  dis¬ 
tances  the  body  appears  as  a  point  source,  the  velocity  potential  according 

22  22  -1/2 

to  linear  inviscid-flow  theory  is  proportional  to  €  (x  -B  r  )  for  x  >  Br. 
Therefore  u  =  as  x(Br)  *  -  1  0.  A  solution  of  the  form 

(2.44),  which  includes  the  effects  of  viscosity,  can  match  as  t  oo  with  this 
potential-flow  solution  if  a  =  3/4;  this  is  equivalent  to  the  results  of  Refs. 

5  and  6.  The  solution  for  g(t)  can  be  written  in  terms  of  confluent  hypergeo- 
metric  functions,  whose  properties  are  known^^.  The  result  is  exponentially 


Abramowitz,  M.  and  Stegun,  I.  A.  (eds* )»  Handbook  of  Mathematical 
Functions.  National  Bureau  of  Standards,  1964. 


small  as  t  -  oo  (upstream)  and  is  of  order  t  as  t  +co  (downstream). 
The  integral  of  f(t)  over  the  range  -  oo  <  t  <  oo  is  found,  by  direct  integra¬ 
tion  of  the  differential  equation  (2. 45),  to  be  zero.  The  second  term  in  u  is 

t 

the  simple  dipole  solution  for  or  =  1,  which  is  exponentially  small  both  up¬ 
stream  and  downstream;  the  nonlinear  term  is  still  of  higher  order.  If 

2 

there  were  no  net  mass  outflow  of  order  e  from  the  neighborhood  of  the 
body,  the  solution  for  a  =  1  would  be  the  leading  term  in  u. 

Transonic  Flows 

The  definition  (2.11)  of  rj  and  the  expansion  (2.12)  for  u  could  have  been 

modified  so  as  to  show  explicitly  the  dependence  on  the  parameter  B.  The 

orders  of  magnitude  for  r  and  u  would  again  be  found  by  specifying  that  a 

nonlinear  term  must  appear  in  the  differential  equation  (2.14)  and  that  the 

solution  for  u  must  match  asymptotically  with  the  solution  (2. 10).  If 

is  then  allowed  to  decrease  toward  one,  so  that  B  -*■  0,  the  nonlinearity 

-4  3  4  -2 

appears  for  r  =  0(e  B  ),  with  u  =  0(e  B  ),  whereas  the  linearized  flow 
description,  in  terms  of  a  source  distribution  along  the  axis,  gives  the  solu¬ 
tion  for  r  =  0(B  V  These  distances  are  no  longer  distinct  when  B  -  O(e), 
corresponding  to  the  usual  transonic  small -disturbance  theory  derived  with 
the  similarity  parameter  K  =  (M^-D/e^  held  fixed  as  €  0  and  M  —  1. 

For  B  =  0(e),  then,  u  =  O(e^)  and  the  first  approximation  satisfies  the  non-¬ 
linear  transonic  small-riisturbance  equation  in  the  variables  x  and  r  =  er. 
The  ratio  of  second-order  to  first-order  terms  is  O(A^),  where  =  (e^Re) 
for  the  transonic  limit;  the  shock-wave  thickness  is  also  of  this  order. 

At  small  values  of  r  such  that  r  =  0(A  ),  the  shock  wave  is  no  longer  one- 
dimensional  but  is  described  instead  by  the  "viscous-transonic"  equation, 
which  is  the  transonic  small -disturbance  equation  augmented  by  the  diffusion 

term  proportional  to  u  .  At  large  values  of  r,  again  as  B  —  0,  the  previous 

XX  -3/4  1/4 

solutions  for  shock-wave  location  and  thickness  give  |  =  0(eB  r  ) 

-1  -1/4  -1  3/4  Si 

and  A  =  0( e  B  Re  r  ).  These  two  lengths  are  no  longer  distinct  if 

r  *  0(e  o*  Re  j  and  4  =  0(e  B  Re  /  In  the  fully  viscous  region  defined 

in  this  way,  the  velocity  is  u  *  0(e  ^B  Re  For  B  =  O(e),  the  solution 
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for  r  =  0(1)  no  longer  has  the  form  shown  by  Eqs.  (2.15)  through  (2.17),  but 

for  r  »  1  the  quantities  %  and  A  have  the  same  power-law  behavior  as 

31  S 

for  B  >>  €,  now  with  unknown  numerical  factors.  The  fully  viscous  region 
for  B  =  O(c)  is  then  defined  by  r  =  Q(A*)  and  £  =  0(A^2),  and  can  again 
be  described  by  the  modified  form  of  Burgers'  equation  required  for  axi- 
symmetric  flow. 

If  the  free-stream  speed  is  exactly  sonic,  the  inviscid-flow  solution 

as  r  -*  oo  has  the  self-similar  form  u  =  e^r "^^ffx/r^^).  Then 

uu  =  0(e^?  an<i  Re  =  0(Re  *c2r  as  r  ■*  oo.  There  now  is  a 

x  xx  -2Jz  /l 

single  shock  wave  having  thickness  A  =  0(Re  c  r  ),  whereas  the 

j4/7 

inviscid-flow  solution  implies  x  =  0(r  ).  By  either  comparison  it  is  seen 

that  the  fully  viscous  region,  described  by  the  viscous-transonic  equation, 

~  -7  /2 

appears  for  r  =  0(A  ).  The  self --similar  solution  then  is  correct  for 

~  -7/2  1  2 

1  «  r  «  A^  .  If  now  B  is  increased  from  zero,  a  linear  term  B  u 
t  x 

will  first  contribute  to  the  description  of  the  fully  viscous  region  when 
B2  is  no  longer  small  in  comparison  with  u;  for  this  condition 

B/e  =  0(A^2).  The  previous  result  r  =  0(e^B  *Re2)  for  the  distance 

*  3/2 

to  the  fully  viscous  region  is  still  correct  for  B/e  as  small  as  0(A  ), 

^  -7/2  3/2  * 

but  is  replaced  by  r  =  0(A  )  for  B/e  «  A  .  The  orders  of  mag- 

t  c 

nitude  for  the  locations  of  regions  in  which  different  approximations  are 

required  are  shown  schematically  in  terms  of  Mach  number  in  Fig.  3. 

For  subsonic  flow  with  M  slightly  smaller  than  one,  the  simi- 

2  °°  2 

larity  parameter  is  K  =  (1-M  )/e  ,  and  the  transonic  small -disturbance 

°°  ^  2  2 

equation  is  obtained  in  coordinates  x  and  r  if  (l-M^)/e  =  0(1).  For 

r  »  1,  the  flow  is  purely  subsonic,  and  the  body  (if  closed)  appears  as 
-2  -  3 

a  doublet  with  e  r  u  equal  to  a  function  of  x/?.  The  linear  term 
2  2  2 _ 4 

(1-M  )u  is  then  of  order  (1-M  )e  r  ,  whereas  the  diffusion  term  is 
proportional  to  Re  *  =  0(Re-^e^r .  That  is,  the  diffusion  term 

decays  faster  and  never  appears  in  the  first  approximation  as  r  -♦  00. 

If,  instead,  0  <  1  -  M^  «  e2,  the  first  approximation  for  r  =  0(1)  is 


the  same  as  for  M  =1,  and  (1-M^  ^  =  0((1-M^  )e^r  as  r  oo. 

oo  oo  x  oo 

This  term  is  no  longer  small  in  comparison  with  the  nonlinear  term 

2  ^  -7  /6 

when  u  and  (1-M  )  are  of  the  same  order,  for  r  =  0(K  /  ).  The 

00 

shock  wave  then  extends  to  a  distance  r  =  0(K  )  and  its  thickness 

-4/7 

remains  small  in  comparison  with  Ax  =  0(r  )  provided  that 

K  ^ ^  «  €^Re^^,  i.e.,  (1-M^  ) / »  A  \  In  other  words,  the 

®  2  *  2  -3 

doublet  describes  the  far  field  for  (1-M  )/e  »  A  and  a  fully  viscous 

oo  t 

cw 

region  appears  for  large  values  of  r  only  in  the  very  narrow  Mach 

2  -4  -3 

number  range  (1-M^)  =  0(e  Re  ). 

Sample  numerical  values  indicate  that  viscous  effects  are  very  small, 
as  expected,  and  in  some  respects  become  less  important  as  the  Mach  num¬ 
ber  decreases  toward  one.  A  convenient  pair  of  values  is  €  *  0.1,  Re  =  10^. 

4 

At  supersonic  speeds,  nonlinear  effects  become  essential  for  rj  =  e  r  =  0(1), 

4 

i.e.  ,  for  r  =  0(10  ),  At  this  distance  the  shock-wave  thickness  is  O(A), 

4-1  1 

where  A  -  (€  Re)  =  0.01;  the  fully  viscous  region  appears  for  rj  =  0(A"^), 

8 

i.e*,  for  r  =  0(10  ).  At  M^  =  1 ,  the  shock -wave  thickness  for  r  =  6r  s  0(1) 

2  - 1  -4 

is  0(Afc),  where  A.  =  (e^Re)  =  10  ,  and  the  fully  viscous  region  appears 

7 / 2  t  15 

for  *r  =  0(A^  ),  i.e.  ,  for  r  =  0(10  ).  The  solutions  for  M^  =  1  are  correct 

in  the  extremely  narrow  Mach-number  range  | 0(A  /*),  i.  e.  ,  |M^  -  1  | 
14  t  oo 

=  0(10  ).  Some  of  these  features  are  apparent  in  Fig.  3. 
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Figure  3.  Mach-number  ranges  for  various  approximations. 


3.  LOCATION  OF  A  SHOCK  WAVE  IN  LOW  REYNOLDS  NUMBER  FLCW 


T.  C.  Adamson,  Jr. 


Introduction 

It  is  well  known  that  in  numerical  solutions  of  supersonic  flow  fields 
the  location  of  a  shock  wave  is  known  only  to  within  three  or  four  mesh 
points  when  the  so-called  shock  capturing  technique  is  used.  That  is,  the 
shock  wave  appears  not  as  a  discontinuity,  but  as  a  thick  region  in  which 
flow  properties  vary  relatively  rapidly.  Often  this  is  not  a  serious  drawback 
because  only  a  single  shock  wave  occurs  and  it  is  not  necessary  to  locate 
the  shock  wave  in  the  flow  field  more  accurately  in  order  to  calculate  forces 
on  the  body  to  the  desired  accuracy.  The  fact  that  the  pressure  distribution 
at  the  intersections  of  the  shock  wave  and  the  wall  is  spread  out  over  a 
region  large  compared  to  the  thickness  of  the  shock  wave  is  not  in  itself 
incorrect  in  view  of  the  fact  that  this  is  what  occurs  in  the  interaction  re¬ 
gion  caused  by  the  intersection  of  a  shock  wave  and  a  boundary  layer;  if 
the  extent  of  this  expanded  pressure  distribution  is  not  too  much  different 
from  that  of  the  interaction  region,  the  errors  incurred  may  not  be  prohib¬ 
itive.  However,  before  accurate  results  using  known  solutions  in  the  inter¬ 
action  region  can  be  calculated,  it  is  necessary  to  locate  the  shock  wave 
accurately.  Moreover,  in  flows  over  complex  bodies  and  in  internal  flows, 
where  complicated  shock  wave  patterns  and  reflections  are  important,  such 
a  lack  of  preciseness  in  location  of  the  shock  wave  does  not  suffice. 

One  of  the  methods  for  calculating  the  location  of  the  shock  wave 
accurately  is  that  of  shock  fitting  (e.g.,  see  reference  1).  Here,  at  each 
iteration  of  the  solution,  the  shock  wave  is  positioned  such  that  the  proper 
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jump  conditions  are  met  and  the  solution  on  either  side  of  the  shock  wave 
merges  with  the  known  numerical  solution.  Unfortunately,  the  logic  in¬ 
volved  in  this  calculation  is  complicated  for  even  a  single  shock  wave  and, 
at  least  in  the  present  form,  too  unwieldy  to  consider  for  complex  shock 
wave  patterns.  It  is  clear  that  a  simple  method  for  locating  the  shock  wave 
would  be  of  considerable  value. 

This  section  covers  the  initial  investigation  of  a  method  different  from 
those  considered  before,  for  locating  a  shock  wave.  It  is  applicable  in 
those  cases  where  mass,  momentum,  and  energy  are  conserved  in  the 
"thick"  shock  wave  given  by  the  numerical  solutions.  The  essential  idea  is 
that  the  thick  shock  corresponds  to  a  shock  wave  in  a  relatively  low  Reynolds 
number  flow.  An  asymptotic  solution  for,  say,  the  distribution  of  velocity 
through  a  shock  wave  is  to  be  compared  with  the  numerical  solution  to  give 
the  values  of  the  necessary  parameters.  Then,  the  analytical  solution  is  to 
be  used  to  locate  the  shock  wave  in  the  actual  relatively  high  Reynolds  num¬ 
ber  flow  by  simply  using  the  solution  valid  in  the  limit  as  Reynolds  number 
tends  to  infinity.  In  actuality,  because  the  location  of  the  sonic  line  within 
the  shock  wave  corresponds  to  the  location  of  the  shock  wave  in  the  real  flow 
case,  where  the  shock  is  very  thin,  the  calculation  reduces  to  finding  the 
dependence  of  the  location  of  the  sonic  line  on  the  Reynolds  number.  As  it 

V 

turns  out,  in  the  flow  problem  considered  there  is  a  range  of  parameters 
for  which  the  sonic  line  (within  the  shock  wave)  is  independent  of  the  Rey¬ 
nolds  number  to  the  accuracy  desired.  In  this  event,  the  shock  wave  is 
very  easily  located  by  locating  the  sonic  line.  For  smaller  Reynolds  num¬ 
bers,  the  equation  derived  for  the  first  order  correction  to  the  location  of 
the  shock  wave  indicates  an  apparent  dependence  upon  Mach  number.  How¬ 
ever,  this  equation  contains  a  term  which  has  not  been  evaluated;  there  is 
a  possibility  that  the  final  equation  will  be  independent  of  Reynolds  number. 

In  order  to  test  the  utility  of  the  proposed  method,  it  is  necessary  to 
choose  a  flow  problem  for  which  solutions  in  the  form  of  numerical  compu¬ 
tations  can  be  found  relatively  easily  and  for  which  analytical  solutions  are 
known.  The  test  then  consists  of  comparing  the  shock  locations  found  using 
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on  the  one  hand  the  numerical  solution  in  conjunction  with  the  shock  wave 
locating  method  and  on  the  other  hand  the  analytical  solution.  The  test  prob¬ 
lem  chosen  is  that  of  steady,  transonic  flow  through  a  2 -dimensional  channel 
arbitrary  wall  shape.  Solutions  for  inviscid  flow  upstream  and  down¬ 
stream  of  the  shock  wave  are  known  in  the  form  of  asymptotic  expansions 
uniformly  valid  to  the  second  order,  1  f  ’  both  for  steady  and  unsteady 
flow.  It  is  necessary  to  modify  these  solutions  somewhat  and  to  add  a  higher 
order  term  and  solutions  for  the  structure  of  a  shock  wave  to  complete  the 
relations  needed  here.  Since  only  the  flow  downstream  of  a  sonic  nozzle  is 
really  needed  for  testing  purposes,  the  numerical  computations  should  not 
be  very  difficult. 

In  this  report,  then,  we  first  present  the  analytical  solutions  for  the 
inviscid  flow,  to  be  used  as  a  basis  for  comparison.  In  these  solutions,  the 
shock  wave  appears,  of  course,  as  a  discontinuity.  Next,  we  consider  the 
same  channel  flow,  but  now  for  the  case  where  the  Reynolds  number  is  rela¬ 
tively  low,  but  still  large  compared  to  unity.  In  this  case,  the  shock  is  rela¬ 
tively  thick,  and  it  is  this  solution  which  is  to  be  compared  with  the  numerical 
solution  to  find  the  equivalent  Reynolds  number,  and  other  necessary  parame¬ 
ters.  It  is  also  this  solution,  then,  which  is  to  be  used  to  find  the  change  in 
location  of  the  shock  wave  as  Reynolds  number  becomes  very  large;  the 
change  is  used  in  correcting  the  numerical  solution. 

It  should  be  noted  that  the  method  of  computation  used  has  much  to  do 
with  the  form  of  ,fthick  shock  wave”  found  in  numerical  computations.  That 
is,  with  a  straightforward  finite  difference  scheme,  one  might  find  a  monotonic 
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decrease  in  velocity,  for  example,  which  resembles  the  distribution  found  in 
a  shock  wave;  with  other  formulations,  an  oscillation  may  be  found  in  the 
velocity  in  the  region  of  the  shock*  In  the  latter  case,  the  oscillations  can  be 
removed  with  the  addition  of  an  artificial  viscosity  and  it  is  supposed  that 
this  would  be  done  if  the  present  method  for  locating  the  shock  wave  were  to 
be  used. 

Problem  Descriptions  -  Solutions  for  In  viscid  Flow 

We  consider  a  steady  transonic  flow  in  a  two-dimensional  channel.  A 
sketch  illustrating  the  coordinate  system  and  notation  used  is  shown  in  figure 
(1),  The  specific  case  chosen  here  is  that  where  the  flow  accelerates  through 
a  sonic  throat  to  a  supersonic  velocity;  at  some  point  downstream  of  the  throat, 
the  channel  walls  become  parallel  to  each  other  and  terminate  in  a  plenum 
chamber,  at  which  a  back  pressure  is  applied  such  that  a  shock  wave  forms 
in  the  channel. 

The  gas  is  assumed  to  follow  the  perfect  gas  law  and  to  have  constant 
specific  heats.  The  flow  upstream  of  the  shock  wave  is  isentropic,  and  the 
shock  itself  is  weak  enough  that  a  velocity  potential  may  be  used  to  the  order 
desired.  Coordinates  x  and  y  are  made  dimensionless  with  respect  to  the 
throat  half  width,  L,  and  velocity  components  with  respect  to  the  undisturbed 
critical  sound  speed,  a*  .  (Overbars  denote  dimensional  quantities.)  The 
pressure,  p,  density,  p,  and  temperature,  T,  are  made  dimensionless  with 
respect  to  their  critical  values,  and  the  enthalpy,  H j  is  referred  to  a*^. 

The  wall  shape  is  written  as  follows,  for  symmetric  channels,  where 
€  «  Is 

yw  *  +  (1  +  e2f(x))  (3.  i) 

where  f(0)  =  f?  (0)  =0.  Thus,  may  be  related  to  the  radius  of  curvature 
of  the  walls  at  the  throat,  x  =  0. 

Composite  solutions,  uniformly  valid  to  order  may  be  written  for 
the  flow  under  consideration  as  follows:^’  ^ 
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2 

1  +  e  u^  +  e  u^  +  . . . 

(3.2a) 

2  ,  5/2  * 

e  v,  +  e  5  +  •  •  • 

2  y 

(3.2b) 

1  -  eYUj  -  €  2  Y  u2  +  ... 

(3.2c) 

,  2.  .  (Y-l)  2,  , 

1  "  €uj  _  €  ^u2  +  ’  u.^)  +  . .  . 

(3. 2d) 

u2 

1  -  c(Y-l)  -  e2(\-l)(u2  +  y )  +  ... 

(3.  2e) 

where  v  =  c  /c  is  the  ratio  of  specific  heats  and  where 
P  v 


U1  =  -  /y+D  £(X)  +  Cw 


U2  =  2”  +  +  =  r  y 

x 


h  =  -  f  [f  +  u?(2\  -3)]  +  — 
x  6  L  1  J  Uj 


*  •  ..,>/*  s  uu 


5  *  «o[,Y+Ilai0| 


cp  ,  .,n 

—  cos  (n  7T  y)  • 


n=  1  (tt  n) 


exp  {  -nir  x  /[(y+Du^q]1/2}  x  >  x{ 


?  =  0 


x  <  x 


*  /  v  -1/2 

X  =  (x  -  X  )  € 

SO 


(3.3a) 

(3.3b,  c) 

(3.  3d) 


(3.  3e) 

(3.  3f ) 

(3.  3g) 


In  eqn.  (3.  3a),  is  an  arbitrary  constant  determined  by  the  value  of  the 
velocity  at  the  throat;  i.  e.  if, there,  the  flow  is  subsonic  or  supersonic, 

*  0.  The  constant  C  is  arbitrary  and  may  have  different  values  up¬ 
stream  and  downstream  of  the  shock  wave.  In  the  equations  above, 
f  =  df/dx,  etc.  Also,  f'^  and  u^  are  f  and  the  positive  value  of  u^. 
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respectively,  evaluated  at  x  ,  the  zero -order  approximation  to  the  location 

s  o 

of  the  shock  wave,  x  ,  where  x  is  expanded  as  follows: 

s  s 


*3  =  *30  +  e*,l  +  eV2*.3/2(/)  +  '2*s2W  +  ••• 


(3.4> 


That  is,  dx  /dy  =  [v]/[u]  =  0(e^^/e)  =  where  the  square  brackets 

s 

indicate  the  jump  in  the  quantity  enclosed  across  the  shock  wave  at  the  point 

(2) 

y  and  x  (y),  but  it  can  be  shown  that  there  is  a  term  of  O(c)  in  the  expansion. 

Hence,  in  eqn.  (3.4),  x  and  x  .  are  constants  for  steady  flow.  It  should  be 

s  o  s  i 

noted  that  in  this  inviscid  flow  solution,  the  shock  wave  is  a  discontinuity. 

* 

The  quantity  ?  which  appears  in  eqns.  (3.2b)  and  (3.  3b)  and  which  is 

defined  in  eqns.  (3.  3e)  and  (3.  3f )  is  an  additional  potential  function  needed  in 

(2) 

the  neighborhood  of  the  shock  wave.  That  is,  it  can  be  shown  that  the  chan¬ 
nel  flow  solutions  satisfy  the  jump  conditions  across  the  shock  wave  to  first 

order  but  not  in  second  and  higher  order  terms.  Hence,  it  is  necessary  to 

1/2 

consider  a  region  of  order  e  in  thickness  in  the  neighborhood  of  the  shock 

wave;  in  this  region  the  flow  must  adjust  from  the  conditions  required  by  the 

shock  wave  to  those  associated  with  the  outer  channel  flow.  In  general,  be- 

* 

cause  the  flow  is  supersonic  upstream  of  the  shock  wave,  ?  =0  there  and 

an  adjustment  region  is  needed  only  downstream  of  the  shock  wave.  There, 

*  * 
x  (defined  in  eqn.  (3.  3g)  is  0(1)  and  as  shown  in  eqn.  (3.  3e),  ?  goes  to  zero 

$ 

exponentially  as  x  increases.  As  mentioned  previously,  the  solutions  shown 
in  equations  (3.3)  are  composite  solutions;  the  solutions  in  the  inner  region 
have  been  included  and  the  solutions,  as  written,  are  uniformly  valid  to  order 
throughout  the  channel. 

The  location  of  the  shock  wave  may  be  found  in  either  of  two  ways.  On 

the  one  hand,  the  composite  solutions  may  be  used  in  conjunction  with  the  jump 

conditions  across  the  shock  wave.  On  the  other  hand,  conservation  of  mass, 

y 

*  w 

enforced  by  equating  /  p  udy  evaluated  upstream  and  downstream  of  the 
shock  wave,  may  be  used.  The  former  method  involves  the  evaluation  of 
higher  order  solutions  in  the  inner  region  if  higher  order  approximations  to 
the  location  of  the  shock  wave  are  desired,  and  so  the  latter  method  is  employed 


here.  As  pointed  out  in  reference  2,  it  is  necessary  to  calculate  the 

4 

jump  in  entropy  across  the  shock  wave  and  pu  to  order  e  if  one  wishes 

to  calculate  x  _.  In  addition,  u^#  the  third  order  outer  solution  for  u 
si  3 

must  be  known;  again,  if  the  mass  flow  is  evaluated  far  enough  upstream 
and  downstream  of  the  shock  wave,  then  inner  solutions  are  unnecessary. 
The  solution  for  u^  shown  below  has  been  found  using  the  methods  des¬ 
cribed  in  reference  (2).  The  required  equations  for  pu  and  the  jump  in 
entropy  across  the  shock  wave,  also  shown  below,  are  essentially  those 
given  in  reference  2;  a  few  corrections  have  been  incorporated: 


u3  *  <V+»CV<24  •  1?  *T  Y*]xx  *  8* 


(3.  5a) 


„  _  (Y+1)f_Z_  u  £n  .  — Li  +  r  -  .  YiilziLi  3 
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=  1  -  €2  u2  -  €3(y+1)[u1u2  +  u3] 
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where  the  entropy,  S,  is  made  dimensionless  with  respect  to  the  gas  constant, 

R. 

If  eqn.  (3.  5d)  is  integrated  at  x  =  0  and  at  x  =  X,  say,  where  X  >  x  , 

2  s 

and  the  results  are  equated,  the  terms  of  order  e  reduce  to  eqn  (3.  3a), 

3  4 

and  the  terms  of  order  e  and  e  give,  respectively. 


"2d 


=  C 


2u 


2l  3 

3  10 


3d 


°3u  “to  -  ^‘./loVo  -  2l,U10C2u 


(3.6a) 

(3.6b) 


where  u,  -  refers  to  the  value  of  u,  (x  )  immediately  upstream  of  the  wave. 
1x0  lx  so 

Thus,  setting  the  constants  C  and  C  ,  equivalent  to  setting  the  pressure 

Ct  U  k  U 

to  second  order  at  points  upstream  and  downstream  of  the  shock  wave,  allows 

one  to  calculate  u1A,  from  eqn.  (3.6a);  therefore,  from  eqn.  (3.  3a),  x 

10  so 

can  be  calculated.  Next,  setting  C  and  C  ,  equivalent  to  setting  the  pres- 

ju  jll 

sure  at  the  same  points  to  third  order,  one  may  calculate  x^  from  eqn. 
(3.6b).  Thus,  with  the  inviscid  flow  solutions  shown  here,  the  position  of  the 
shock  wave  can  be  calculated  to  first  order  accuracy. 


Solutions  Including  Structure  of  the  Shock  Wave 

The  solutions  in  the  previous  section  are  found  under  the  assumption  that 

the  shock  wave  is  a  discontinuity,  and  are  valid  as  long  as  the  thickness  of 

the  shock  wave  is  small  compared  to  the  order  of  the  accuracy  to  which  the 

shock  is  located.  That  is,  if  the  order  of  the  dimensionless  (with  respect  to 

the  throat  half  width),  thickness  of  the  shock  wave  is  given  by  L  ,  then,  as 

indicated  in  the  following  sketch,  x  is  known  to  order  e  as  long  as  L  «  €. 

s  s 


OUt) 


SHOCK  WAVE 
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In  order  to  relate  the  orders  of  L and  e  to  be  used  in  the  present 
application,  it  is  necessary  to  estimate  the  thickness  of  a  shock  wave 
"captured"  in  a  typical  numerical  computation  relative  to  the  character¬ 
istic  length  in  the  problem  under  consideration.  For  example,  in  typi¬ 
cal  computations  of  flows  over  airfoils,  the  shock  thickness  can  vary 
from  roughly  three  percent  (Reference  5)  to  eight  percent  (Reference  6) 
of  the  chord;  in  computations  of  the  flow  through  a  stator  blade,  with 

artificial  viscosity  added,  this  thickness  can  be  as  high  as  12  percent 
(7) 

of  the  chord.  If  the  effects  of  adding  artificial  viscosity  are  neglected 
for  the  present,  it  appears  that  a  good  estimate  of  the  thickness  of  the 
captured  shock  wave  is  five  percent  of  the  characteristic  length  in  the 
flow  direction.  In  the  present  case,  the  throat  width  is  the  character¬ 
istic  length,  i.  e. ,  L  =  0.05(2L),  or 

s 


t  -  _  J.  (2L  _L 

s  L  20  '  L;  10 


(3.7) 


If  the  Reynolds  number,  Re,  (based  on  L  and  critical  flow  conditions) 
and  the  Prandtl  number,  P  ,  are  written  in  terms  of  the  longitudinal 

r  (g) 

viscosity,  then  following  Illingworth 


Yoshihara,  H. ,  Finite  Difference  Procedure  for  Unsteady  Transonic  Flows: 

A  Review,  Transonic  Flow  Problems  in  Turbomachinery,  eds.  T.  C.  Adam¬ 
son,  Jr.  and  M.  F.  Platzer,  Hemisphere  Publishing  Corp. ,  1977,  139  -162. 

^Caughey,  D.  A.  and  Jameson,  A.,  Calculation  of  Transonic  Potential  Flow- 
fields  About  Complex  Three-Dimensional  Configurations,  Transonic  Flow  Prob¬ 
lems  in  Turbomachinery.  eds.  T.  C.  Adamson,  Jr.  and  M,  F.  Platzer, 
Hemisphere  Publishing  Corp.,  1977,  274-293. 

^Farn,  C.  L.  S.  and  Whirlow,  D.  K. ,  Application  of  Time -Dependent  Finite 
Volume  Method  to  Transonic  Flow  in  Large  Turbines,  Transonic  Flow  Prob¬ 
lems  in  Turbomachinery.  eds.  T.  C.  Adamson,  Jr.  and  M.  F.  Platzer, 
Hemisphere  Publishing  Corp.,  1977,  208-227. 

^Illingworth,  C.  R. ,  Shock  Waves,  Modern  Developments  in  Fluid  Mechanics. 
High  Speed  Flow,  Vol.  I.  ed.  L.  Howarth,  Oxford  Press,  1953,  105-137. 
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L  =  - - -  (1  +  'L=-^  )  (3.8a) 

8  Re  e(v+l)  Pr 

*  4k  (3.8b) 

s 

where  equation  (3.8b)  defines  the  parameter  kg,  to  be  used  hereafter  as 

being  of  the  order  of  the  thickness  of  the  shock  wave.  It  is  seen,  from  eqns. 

(3.7)  and  (3.8b),  that  the  value  of  k  corresponding  to  the  estimated  thickness 

s 

of  the  shock  wave  found  in  numerical  studies  is  k  =1/40  and  further,  from 

_  s 

eqn.  (3.8a),  that  for  =  1,  Re  e  a  40  y/(v+l).  Finally,  a  typical  value  for 
e  is  1/10. 

It  appears  from  the  above  estimates,  that  solutions  from  either  of  two 
limit  processes  might  be  useful  in  applications  to  numerical  problems.  The 
two  possibilities  are 

(a)  k  =  0(e2);  Re  =  0(e~3) 

s 

(b)  kg  =  0(0;  Re  =  0(e'2). 

In  case  (a),  setting  k  *  me^  with  m  being  a  constant  of  order  unity,  a  nu- 

s 

merical  value  of  2.  5  for  m  would  result  in  the  desired  relative  values  of  k 

s 

and  €.  In  case  (b),  with  k  =  me,  the  corresponding  value  of  m  is  i/4. 

s 

In  the  following  section,  solutions  for  case  (b),  kg  =  me,  are  con¬ 
sidered.  This  case  was  chosen  for  the  initial  study  because  solutions  are 
valid  for  case  (a)  also,  in  the  present  application.  It  will  be  shown  that  it  is 
possible  to  obtain  the  desired  solutions  in  terms  of  an  integral  which  can  be 
evaluated  using  known  functions. 

(a)  Channel  Flow  Solutions;  k  =  me 

s 

The  solutions  considered  here  are  for  a  flow  in  which  the  longitudinal 
viscosity  is  large  enough  that  its  effects  are  found  in  the  second  order  outer 
channel  flow  solutions.  In  addition,  it  will  be  seen  that  the  shock  wave  is 
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thick  enough  that  the  effects  of  streamline  curvature  are  found  in  second  order 
solutions  for  the  shock  structure;  as  a  result,  the  Rankine -Hugoniot  jump 
conditions  do  not  hold  across  the  shock  wave.  These  two  effects  cause  for¬ 
midable  difficulties  in  obtaining  higher  order  solutions. 

The  governing  equations  are  the  Navier -Stokes,  continuity,  and  energy 
equations  for  steady  flow.  Thus: 


uu  +  vu 
x  y 


(3.9a) 


uv  +  vv 
x  y 


(3.9b) 


uH  +  vH  =  ~  (up  +  vp  )  +  ——  R(l-n 
X  y  P  Y  px  p~e 


(3.9c) 


where 


(5x  u  -  (  7  jjl  -  \il  )  v  ]  +  [|i(u  +  v  )] 

x  3  y  x  y  x  y 


(3.  10a) 


=  lvtvx  +uy)]x  +[Hvy  -  <  JH  -  Ii')ux]y 

=  [  H  ]  +  [  H  ]  +  +  w" ] 

~  xx  ~  y  y  x  y 


(3.  10b) 


-(“|i.-2ji')uv  +(u  +  v  ) 

3  r  x  y  y  x 


(3.  10c) 


The  dimensionless  longitudinal  coefficient  of  viscosity,  Reynolds  number, 
and  Prandtl  number  are  defined  as  follows  where  /!'  is  the  bulk  viscosity  and 
|4  the  usual  viscosity  coefficient; 


sr  4  - 

U  =  JjU  +/U 


(3.  11a) 


£  =  mV(m)*  m  = 


(3.  lib) 
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The  equations  of  state  and  relationship  between  enthalpy*  and  temperature  are, 
respectively,  in  dimensionless  terms. 


P  =  p  T 

(3.  12a) 

T 

H  =  — r 

Y-l 

(3.  12b) 

Finally,  three  equations  which  are  of  importance  in  the  following  analysis 

the  gasdynamic  equations  and  those  governing  the  vorticity  and  entropy! 

are 

2  2  2  2 
(u  -a  )u  +  uv(u  +v  )  +  (v  -  a  )v 

X  y  y 

=  [uR(u)  +  vR(v)  -  (y-l)R(H)] 

(3.  13a) 

p  Re 

uft  +  v  ft  =  -  ft  (u  +  v  )  +  ~  ( T  S  -  T  S  ) 

x  y  xyyxy  yx 

+  j 

^  w  p  X  p  y 

Re 

(3.  13b) 

^(uSx  +  vS  )  =  —  R(H)' 

*  y  Re 

(3.  13c) 

2 

where  a  =  T  and  ft  =  v  -  u  are  the  dimensionless  speed  of  sound  with 

x  y 

respect  to  a*  and  z  component  of  the  vorticity  vector  respectively. 

Following  the  same  methods  employed  in  obtaining  the  in  viscid  flow 

(2 ) 

solutions,  one  can  find  the  outer  channel  flow  solutions  valid  now  for  the 

case  k  =  me.  In  order  to  simplify  the  calculations,  Pr  =  1  and  p.'  =  0 
s 

are  taken  to  be  the  case;  more  general  values  will  add  complexity  but  not 
change  the  fundamental  conclusions.  First,  if  the  stagnation  enthalpy,  H 


is  expanded  as  follows 


U2  +v2 

Ht  =  H  +  “i 


_T_  u  +v 

\-l  2 


(3.  14) 


where  -  f.  ,(0  ^  1»  then  from  the  governing  equations  for  Ht» 


u(h  )  +  v(H.) 

t  x  t  i 


— —  [ul  (u)  +  vR(v)  +R(H)] 
p  Re 


(3.  15) 


one  can  show  that  fi  ^  =  e  '  ,  so  that  T  =  a  can  be  calculated  from  the 
equations: 


(u2  +  v2)  =  +0(€9/2) 


1 3. 16) 


Also,  from  eqn.  (3.  13c),  one  can  show  that  As  =  C(e  )  and  from  eqn.  (3.  13b) 

3  4  3 

that  v  -  u  =  0(k  e  )  =  0(€  )  so  that  to  order  €  a  velocity  potential 

X  y  S 

exists.  Finally,  the  boundary  conditions  are  given  by  the  relation 


vU,  ±yw)  =  ±  u(x,  +yw)  (3.17) 

where,  now,  y^  includes  the  displacement  thickness  of  the  boundary  layer 
for  the  low  Reynolds  number  under  consideration.  Now,  because  an  inter¬ 
action  occurs  at  the  intersection  of  a  shock  wave  and  a  boundary  layer  (as¬ 
sumed  to  be  laminar  here),  it  is  important  to  ascertain  the  effects  of  the 
interaction  on  the  displacement  thickness.  In  order  to  estimate  these  ef¬ 
fects,  the  solutions  for  the  flow  within  the  boundary  layer  in  the  interaction 

(9) 

region,  given  by  Brilliant  and  Adamson  ,  may  be  used  with  the  orders  of 
the  velocity  perturbations  valid  outside  the  boundary  layer  found  in  the 
present  calculations.  It  is  not  difficult  to  show  that  at  least  in  lowest  order. 


(9) 

Brilliant,  H.  M.  and  Adamson,  Jr.,  T.  C. ,  Shock  Wave  Boundary -Layer 
Interactions  in  Laminar  Transonic  Flow,  ALAA  J.  ,  12  (  1974),  323-329. 


the  increase  in  v  due  to  the  effects  of  an  interaction  between  the  shock  wave 

and  the  boundary  layer  is  small  compared  to  the  increase  caused  by  the 

change  in  displacement  thickness  over  the  extent  of  the  interaction  region. 

Moreover,  when  the  solutions  valid  within  the  shock  wave  itself  are  con- 

5/2 

sidered,  it  will  be  seen  that  through  0(e  \  ),  the  expression  for  v  is  sim¬ 
ply  a  continuation  of  the  outer  solution  for  v  and  thus  satisfies  the  undisturbed 

flow  boundary  conditions  at  the  wall  to  the  same  order.  This  is  not  the  case 
3  3 

for  the  term  0(e  );  the  term  0(  €  )  in  the  solution  for  v  within  the  shock  wave 

does  not  satisfy  the  boundary  condition  associated  with  undisturbed  flow  at 

3 

the  wall.  It  is  possible  that  this  could  lead  to  corrections  of  Ofe  )  in  the 
displacement  thickness  downstream  of  the  shock  wave  and  thus  to  corrections 
in  second  and  higher  order  terms  in  the  solutions  for  u  downstream  of  the 
shock  wave;  these  corrections  are  not  considered  here. 

If  eqn.  (3.  16)  is  used  for  in  eqn.  (3.  13a),  eqns.  (3.8)  with  Pr  =  1 
are  used  for  Re  ,  and  the  boundary  conditions  to  each  order  of  approximation 
are  obtained  from  eqn.  (3.  17),  the  following  results  are  found  for  u  and  v 
written,  as  in  the  inviscid  flow  case,  in  terms  of  asymptotic  expansions: 


2  3 

u  =  1  +  eiij  +  e  u^  t  £  ^  +  ... 

2  3 

V  =  €  V.  +  €  V-  +  .  .  . 


(3. 


(3. 


where  u^  is  given  by  eqn.  (3.3a)  and 

...  2 


u. 


+  h 


v2  =  f'y 
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so 


U  f" 

v3  =  {y  +1)  [  ~ 


3 

(y  -y)  +  —  1, 


(3.  19f) 


Thus,  the  solutions  are  similar  to  the  corresponding  solutions  in  the  inviscid 
flow  case,  the  difference  being  the  terms  depending  upon  m.  The  same  nota¬ 
tion  is  used  in  the  two  cases,  for  h  and  g;  there  should  be  no  confusion  since 
either  inviscid  flow  or  viscous  higher  order  flow  solutions  are  employed  ex¬ 
clusively  in  any  given  case.  As  in  the  inviscid  flow  case,  constants  and 
Cg  are  arbitrary  and  may  have  different  values  upstream  and  downstream  of 
the  shock  wave,  denoted  by  subscripts  u  and  d,  respectively.  In  the  integral 
in  eqn.  (3,  19e),  the  lower  limit  may,  of  course,  be  chosen  arbitrarily,  the 
choice  being  reflected  in  the  value  of  C  as  boundary  conditions  are  applied, 

w 

Here,  the  choice  is  x^,  the  position  of  the  shock  wave;  upstream  of  the 

shock  wave,  the  integral  has  a  negative  value  and  downstream  of  the  wave 

a  positive  value.  Since  only  the  lowest  order  value  of  the  integral  is  required 

here,  x  is  replaced  with  x  .  Again,  in  eqn,  (3.  19e)  k  refers  to  the  expo- 
s  so 

nent  of  the  temperature  in  the  expression  for  the  viscosity,  taken  to  be 

£T  =  (3.20) 
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It  may  be  noted  that  if  k  =  0(e  ),  then  u  would  be  independent  of  m 

s  b 

and  only  linear  terms  in  m  would  occur  in  uQ. 

^  1/2 

In  order  to  consider  the  flow  in  the  region  of  order  e  in  thickness 

downstream  of  the  shock  wave,  it  is  necessary  to  know  the  flow  variables  just 

downstream  of  the  wave*  In  the  inviscid  flow  case,  the  jump  conditions 

across  the  shockwave  supplied  these  conditions.  In  the  present  case,  it  is 

not  clear  that  the  jump  conditions  hold;  their  validity  must  be  checked.  In 

this  regard,  the  continuity  equation  (pu)  +  (p  v)  =0  provides  a  simple 

x  y 

test.  Thus,  if  this  equation  is  integrated  across  the  shock  wave,  then  the 

} 

jump  in  p  u,  indicated  by  [pu],  is 

[pu]  =  0(ks(pv)y) 

2  2  3 

since  p  =  0(1)  and  v  =  0(e  ),  [pu]  =  0(k  €  )  =  0(e  ).  Thus,  the  divergence 

s 

of  the  streamlines  within  the  shock  wave  begins  to  affect  [p  u]  in  third  order 
terms.  However,  from  eqn.  (3.  5d),  it  is  seen  that  the  third  order  terms 
in  p  u  involve  second  order  terms  in  u.  Hence,  for  kg  =  O(e),  the  jump 
conditions  break  down  in  second  order  terms;  the  shock  wave  is  a  Hugoniot 
shock  to  first  order  only.  As  a  result,  solutions  valid  within  the  shock  wave 
itself  must  be  used  to  continue  the  solutions  across  the  shock. 


(b)  Solutions  in  Inner  Structure  Region 

and  Calculations  of  x.  and  x,  , 
_ io _ t  1 

The  following  sketch  shows  the  various  regions  under  consideration. 


and  the  orders  of  thickness  of  the  inner  regions.  Thus  the  outer  regions  have 

lengths  of  0(1).  The  inner  adjustment  region,  employed  also  in  the  inviscid 

flow  calculations  mentioned  previously,  has  a  length  scale  of  order 

the  inner  structure  region  has  a  thickness  of  order  k  in  general,  so  that  for 

s 

the  present  calculation,  its  thickness  is  0(e). 

Because  the  length  scale  in  the  inner  structure  region  is  the  order  of 

the  thickness  of  the  shock  wave,  the  independent  variables  in  this  inner  region 
+  + 

are  taken  to  be  x  and  y  ,  where 


+  x-Xj  (y) 


+ 

y  -  y 


(3.21a,  b) 
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and  where  x^  (y)  is  the  location  of  the  sonic  line  within  the  shock  wave.  That 

is9  in  general,  the  sonic  line  is  chosen  to  represent  the  position  of  a  shock 

wave  with  finite  thickness,  and  the  distributions  of  velocity,  pressure,  etc. 

are  written  with  respect  to  the  sonic  line.  In  the  limit  as  Re  oo  and 

k  0,  it  is  clear  that  x.  =  x  ,  the  location  of  the  shock  wave  in  the  inviscid 
s  is 

flow  solution.  Hence  it  is  seen  that  a  similar  expansion  may  be  used  for  ; 


i .  e . , 


+  ex(,  +  c  3//2j 


xi  =xio  +  €Xil  +  C  Xi  3/2(y)  +  €  X12(y)  +  (3* 

The  goal,  then,  is  to  find  x^  ^  and  in  terms  of  m;  as  m  —  0 

(k  /e  —  0)  x.  —  x  and  x,  -►  x  .. 
s  i  o  so  i  1  si 

If  the  outer  solutions  upstream  and  downstream  of  the  shock  wave  are 
expanded  about  x^  and  written  in  terms  of  inner  variables,  the  following  re¬ 
sults  are  found: 

u  =  i  +  eu^x^  Q)  +  e2  {(x^  L  +  mx+)ulx(xi  q)  +  u2(x^  0' 

+  €  5/2{xi  3/2Ulx(xi  o>  +  (xi  1  +  ^  C  x*  (0'  y)} 

,  ^  V 

3  (x^+mx  )  + 

+  e  {x  a  (x  )  + - - -  u  (x  )  +  (x  +  mx  )u  (x  ,  y) 

12  lx  io  2  lxx  to  11  2x  1  o 


(3.22) 


+  *13/2  W0'*  * 


(x  titix4)  „ 

— - - 


+  u3(xi  Q,  y)  +  77  (0,  y)}  + 


(3.23a) 


5/2  * 


v  =  e  f  (x.  )y  +  e  ?  (0,  y) 


+  e3{(x  +mx+)(f"(x  )+?  .  (0,  y))  +  v  (x  ,y)}+...  (3.23b) 

II  I  o  yx*  o  i  o 

where  the  expansion  for  x^  (eqn.  3.  22)  has  been  used  and  where  it  has  been 
anticipated  that  a  composite  solution  similar  to  that  used  in  the  inviscid  flow 
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problem  will  be  necessary.  Thus,  ?  and  r]  are  the  second  and  third  order 

potential  functions,  respectively,  which  are  needed  in  the  inner  adjustment 

region  downstream  of  the  shock  wave;  they  are  both  zero  upstream  of  the 

wave.  Anticipating,  further,  that  the  lowest  order  terms  in  the  relations 

for  the  locations  of  the  sonic  line  and  the  inviscid  flow  shock  wave  are  the 

same,  x  =  x  ,  it  may  be  noted  that  the  notation  introduced  previously 
so  x  o 

for  functions  of  x  will  be  used  in  evaluating  equations  (3.  2  3)  in  what  follows; 
s  o 

i.  e. ,  f(x^  q)  =  fQ,  u^(x^  )  =  upstream  of  the  wave,  etc. 

In  view  of  eqns.  (3.23)  and  (3.21),  it  is  seen  that  in  the  inner  structure 
region,  the  expansions  for  u  and  v  and  the  derivative  transformations  are 


,  .  t  c  T 

u  =  l  +  cu1  +eu^  +  € 


'“5/2  + 


2  +  5/2  +  .  3 

=  £  V2  +  £  v5/2  +  £  v: 


3  x 


k 

s 


3_ 

3  y 


where 


3/2  , 

€  x' 


t  3/2 


e2x< 


i  2 


(3.24a) 

(3.24b) 


(3.24c,  d) 


(3.25) 


If  eqns.  (3.24)  and  an  expansion  forHt  are  substituted  into  eqn.  (3.15),  with 

— * 

Pr  =  1  and  ju '  =0,  it  can  be  shown  that  within  the  shock  wave 


H 


t 


T 
V  -I 


+ 


+  v2) 


Jx+iL 

2(y  -1) 


+  0(c 


(3.26) 


so  H.  =  constant  to  the  order  retained;  thus,  expansions  can  be  obtained 
1  2 

for  T  s  a  in  terms  of  the  velocity  components.  In  the  same  manner,  eqns. 

(3.13c)  and  then  (3.13b)  can  be  used  to  show  that  within  the  shock  wave,  for 

the  case  considered  (k  =  O(e),  Pr  =  1,  =  0) 

s 

AS  =  0(e2)  (3.27a) 

n  =  v  -u  =  0(€5//2)  (3.27b) 

x  y 
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5/2 

It  may  be  noted  that  £2  is  at  most  0(e  );  i.  e.  this  is  an  order  estimate  of 

5/2 

£2  and  not  a  solution  so  that  the  term  of  order  c  could  be  zero  in  which 
case  £2  would  be  of  even  higher  order.  If  eqns.  (3.24)  are  substituted  into 
eqn.  (3.27b),  it  is  found  that 


+  + 

v  +  =  mu  + 

2  x  ly 


5/2*+ 


(3.28a) 


(3.28b) 


v  =  mu  -  x'  u 
3x+  2y+  iZlx+ 

If,  as  will  be  shown  later,  u,  =  u*(x+),  then  v*  =  v+(y+);  also,  the  equation 

+  Z  Z  + 

for  v^2  *S  lntegrated.  If  these  two  results  are  written  for  x  —>00 

and  x  -  00  and  matched  with  the  corresponding  outer  solutions  calculated 

from  eqn.  (3.2  3b),  it  is  easily  shown  that 


T  + 

v  =  r  v 
2  o7 


(3.29a) 


tf  5/2  =  ’ti3/2(u10‘Ul) 


i  3/2 


(3.29b) 

(3.29c) 


It  is  seen  that  eqn.  (3.29c)  states  that  to  lowest  order,  dx^  /dy  =  [v]/[u] 
where  the  brackets  indicate  the  jump  in  the  quantity  enclosed  across  the  shock 
wave  at  the  point  in  question.  Of  course,  the  shape  of  a  discontinuous  shock 
wave  is  given  by  the  same  equation.  Thus,  even  for  this  thick  shock  wave 
case,  the  differential  equation  (lowest  order)  for  the  location  of  the  sonic 
line  within  the  shock  wave  is  the  same  as  that  for  the  location  of  a  discontin¬ 
uous  shock  wave  (i.  e. ,  shock  in  the  limit  as  Re  —  00),  with  the  same  condi¬ 
tions  upstream  and  downstream  of  the  wave.  This  is  indicated  by  the  fact 
that  x^  2/2*  ecln*  (3.29c),  is  independent  of  m.  Of  course,  it  still 
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remains  to  calculate  the  value  of  x^  at  y  =  0,  i.  e.  x^Q/x^  ^  and  x^  3/2^ 
before  one  can  say  that  the  location  of  the  sonic  line  gives  the  location  of  the 
shock  wave  independent  of  Re. 

If  the  expansions  for  u  and  v,  eqns.  (3.24a,  b)  and  (3.29a,  b)  and  the 
inner  variable  transformations,  eqns.  (3.24c,  d)  are  substituted  first  into 
eqn.  (3.26)  to  obtain  a  ,  and  then  into  eqn.  (3.  13a),  the  following  governing 
equations  are  found  for  the 


U1  u  +  S  a  +  + 

lx  lx  x 


(3.  30a) 


/  +  +  U1  \  H^O  +  (**  +  \  +  + 


(3.  30b) 


(ul+u5 


5/2  K  '  h  +  1)  *(  3/2  ("l0  ‘  U1  *  °  “  5/2;t+x+ 


/  ,  +> 2  .  +  2 
(  +  +  .  ,U2  >  .  (U1  1  + 
K  U3  +  ~2  +  T“  “2  . 


1  7  +  m  l  + 

2  U2  )  +  '  (y+1)  \  3  + 

x  y 


(3.  30c) 


(x'  ) 

'  r  + 


\  *  •  *1  -1)ultr  o)  *  %  +  ♦  +  fcl  “+  +) 


2x  'x 


+  \  +  +  +  + 

+  ("2“,  t)  +  +2u 1  “  ♦  ++2a2U,  +  + 

v  lx  7x  2x  x  lx  x 

(v  +1)  +2  +  „+■  + 

+  2_ul  U,  +  +  •S2a,  +  +  t2ul  K'\  ♦  )  ♦ 

lx  X  lx  X  v  lx  7x 


+  7  + 

+  “l  a,  +  + 

lx  X 


(3.  30d) 
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where  is  the  lowest  order  term  in  the  expansion  for  the  entropy  and  arises 
from  the  expansion  for  p  in  eqn.  (3,13a).  Thus, 


-  1  -AS 
P  =  e 

and  if  eqn.  (3.26)  is  used  to  find  an  expansion  for  T  in  terms  of  the  velocity 
components,  and  if  S  =  S  +  +  . . .  then  it  is  found  that 

O  £ 


Pl  =  -«! 


(3.  31a) 


+  +  (v  -  1)  f  +v2  „  + 

>  =  -  -  11  „  '  (u,  )  +  S. 


(3.  31b) 


Finally,  from  eqn.  (3.13c),  it  can  be  shown  that  the  differential  equation 
for  S  +  is  easily  integrated  to  give 


S2+  .  -  (.  ♦  l)  «  % 


(3.  32) 


where  the  function  of  integration  is  found  to  be  zero  from  matching  consider¬ 


ations. 


The  solution  to  eqn.  (3.  30a)  is  the  well  known  Taylor  solution. 


-  u^  tanh  r 


(3.  33a) 


+  10  +  , 

r  =  -j”  x  ( 

where  the  constant  u_n  in  eqn.  (3.  33a)  results  from  matching  with  the  first 

order  term  in  eqn.  (3.23a)  in  the  limits  x  -*■  -  oo,  x  x,  and  x  —  +  oo, 

<+)  io 
x  x.  .  In  general,  eqn.  (3.  33b)  should  include  a  function  of  integration, 
to 

+  U10  +  ,  r{  , 
r  =  -J—  x  +  G(y). 

However,  if  a  composite  solution  is  formed  from  the  first  order  inner  and 

outer  solutions  and  the  condition  is  enforced  that  at  the  sonic  line 
2  2  2 

u  +  v  =  a  ,  which  reduces  to  u  =  1  for  terms  up  to  first  order,  it  is  easily 
shown  that  G(y)  =  0. 


(3.  33b) 
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The  solution  to  eqn.  (3.  30b)  may  be  written  as  follows: 


.  » ,  m  f'  + 

+  [  2  \2  o  r  .  + 

u.  =  -  -  .  ...  <  —  tanh  r 

2  Wo/  +  1}  2 


(1  +  tanh2  r+) 


/  r*  y, 
' 2  cosh  r  / 


.  ,  ,  ,v  2  In  cosh  r 

+  ktv  -  »»J0  — JT 

cosh  r 


2  4  4 

2  tanh  r  In  cosh  r 

X10V  2  "  .2  + 

cosh  r 


g4(y) 


+  r+  l  *5^ 

tenr  +  — TT  +  “IV 

cosh  r  cosh  r 


Then  eqn.  (3,28c)  may  be  employed  to  find  v ^  .  Thus, 


2mg‘  4  +  +  2m  + 

— r -  r  tanh  r  + -  g'  tanh  r 

»  u_  o 

°U>  10 


(3.  34) 


"'iz")  +S7('r) 


(3,  35) 


where  in  eqns.  (3,  34)  and  (3.  35),  the  g.(y)  are  functions  of  integration, 

4 

Equation  (3,  35)  may,  then,  be  used  to  calculate  v  ,  which  is  needed  in  order 

4*  3y 

to  solve  eqn.  (3.  30d)  for  u^. 

At  this  point,  it  is  convenient  to  use  the  solutions  found  so  far  to  illus¬ 
trate  the  manner  in  which  the  conditions  downstream  of  the  shock  wave  can 
be  found  and  to  calculate  x^  .  That  is,  as  indicated  previously,  the  jump 
conditions  across  a  discontinuous  shock  wave  hold  for  the  present  MthickM 
wave  only  to  first  order.  That  they  do  hold  in  first  order  is  shown  by  com¬ 
paring  the  jump  condition  with  the  first  order  inner  structure  solution,  eqn. 

4 

(3.  33a);  the  inner  solution  indicates  that  as  r  goes  from  -  oo  to  +  oo, 
u^  =  -  u^*  in  agreement  with  the  jump  conditions.  Now,  the  outer 

solutions  evaluated  at  the  shock  wave  position  satisfy  this  first  order  condition 
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because  u^  may  have  either  a  positive  or  a  negative  value.  As  mentioned 
previously,  in  the  case  of  inviscid  flow  the  second  order  jump  conditions 
are  not  satisfied  by  the  outer  solution;  as  a  result  it  became  necessary  to 
consider  the  flow  in  an  inner  adjustment  region.  In  the  present  case,  al¬ 
though  the  Hugoniot  jump  conditions  no  longer  hold  across  the  shock  wave 
to  this  order,  the  second  order  inner  structure  solutions  are  used  to  obtain 
similar  conditions,  which  again  are  not  satisfied  by  the  outer  solutions. 
Hence,  an  inner  adjustment  region  is  still  necessary.  In  summary,  then, 
the  essential  difference  between  the  present  analysis  and  that  employed  in 
the  inviscid  flow  case  is  that  here,  the  higher  order  flow  conditions  down¬ 
stream  of  the  shock  wave,  used  as  the  upstream  boundary  conditions  for 
the  flow  in  the  adjustment  region,  are  found  from  solutions  in  the  inner 

structure  region. 

+  + 

As  r  oo,  then,  u^  becomes: 


+ 
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2 


mfl 
_ o 

(Y  +  1) 


g4(y> 


a. 


10 


fr+-*-oo)  (3.36a) 
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mf’ 

_ o 

(Y  +  1) 


(r+-+oo)  (3.36b) 


If  these  equations  are  matched  with  the  corresponding  outer  solutions,  calcu¬ 
lated  using  eqn.  (3.2  3a),  one  finds  the  following  relations: 


*1  lfl  „ 

uioh H)  '  “io 


-  f  ( 4- 

o  2 


(3.  37a) 


-  *(  lfl  o  +  84^ 
u10 


-  f 


.  (  xi.L) 

O  '  2  6 


10 


"2d 

X10 


-  (0.  y) 
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(3.  37b) 


Hence, 
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(3.  38) 


?**  (0,  y)  =  -  f"Q( y2 
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l,2i  2 

+  — ii 

1  rt  *1  A 


C  -C 
2d  2u 


3  10 


u 


10 


which  is  the  same  result  found  in  reference  2  for  the  inviscid  flow  problem. 

Here,  just  as  in  reference  2,  the  outer  velocity  components  downstream  of 

the  shock  wave  are  expanded  about  the  lowest  order  position  of  the  shock 

*  *  * 

then  £  and  £  are  added  to  these  ex- 

T  y 

V  ' 


wave,  x.  ,  and  thus  in  terms  of  x 
i  o 


pansions.  The  resulting  expressions  are  used  for  the  velocity  components  in 

the  inner  adjustment  regions  in  which  the  flow  variables  change  from  their 

values  immediately  downstream  of  the  shock  wave  to  forms  given  by  the  outer 

*  * 

channel  flow  solutions;  the  additional  potential  function  f  (x  ,  y)  provides 
these  variations.  It  is  not  difficult  to  show  that  in  the  present  problem, 
because  the  expansions  for  the  outer  velocity  components  differ  from  their 
counterparts  in  the  inviscid  flow  problem  only  by  constants,  the  governing 

5$e 

equation,  boundary  conditions,  and  hence  solution  for  £  (x  ,  y)  (eqns. 

* 

(3.  3e)  and  (3.  3f))  are  the  same  as  in  reference  2.  Thus,  since  ?  satisfies 
a  Laplace  equation  and  the  normal  derivatives  on  the  remaining  boundaries 
are  zero,  the  condition 

1 


/  t  *  (0,  y)dy  =  0 


must  be  met.  When  this  condition  is  applied  to  eqn.  (3.  38)  the  results  are 


and 


C  =  C 
2d  2u 


2X  u3 

3  10 


?**(0,y)  =  -r(yZ  -|) 

X 


(3.  39) 


(3.40) 


which  are  the  results  found  in  reference  2.  It  is  worthwhile  to  emphasize 
that  the  results  found  in  eqns.  (3.  39)  and  (3.  40)  are  not  found  if  second  order 
shock  wave  conditions  are  used  to  evaluate  conditions  downstream  of  the  wave; 
again,  the  Hugoniot  jump  conditions  simply  do  not  hold  here. 
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Now,  eqn.  (3.39)  involves  =  u^(x^  an<^  t^lus  maV  used  to  find 
x^  q.  That  is,  by  choosing  the  arbitrary  constants  C^u  and  C^,  one  can 
calculate  u^  and  from  eqn.  (3.  3a),  then,  x^  ,  for  a  given  wall  shape  and 

c  .  Because  eqn.  (3.  39)  is  independent  of  m  and  in  fact  precisely  the  same 

W  (2) 

as  the  equation  found  for  the  infinitesimally  thin  shock  wave  ,  it  is  clear 

that 


x 


i  o 


X 

so 


(3 


i.  e. ,  to  lowest  order  the  location  of  the  sonic  line  within  the  shock  wave 
and  the  location  of  the  shock  wave  under  the  same  flow  conditions  except 
that  Re  —  oo,  are  the  same. 

* 

It  should  be  noted  that  a  composite  solution  involving  ?  and  the  outer 

velocity  components  may  be  formed,  just  as  in  reference  (2);  in  fact,  such 

a  composite  solution  written  to  third  order  in  u  and  thus  involving  a  third 

*  * 

order  potential  function,  rj  (x  ,  y),  has  been  used  to  write  eqn.  (3.23).  In 
this  regard,  eqns.  (3.13c)  and  (3.13b)  can  be  used  to  show  that  Q  =  0(e^'  ^) 

in  the  inner  adjustment  region,  to  confirm  that  a  potential  function  does 
exist  at  least  to  0(e  ). 

In  order  to  find  x^  9  it  is  necessary  to  follow  the  same  procedure  as 

was  used  in  calculating  x-  ,  except  that  now  third  order  rather  than  second 

+ 

order  terms  in  u  are  employed.  That  is,  the  solution  for  not 

volved  in  the  solution  for  x  ;  it  is  necessary,  however,  to  find  u„. 

1  1  +  +  J 
The  governing  equation  for  u  is  found  by  calculating  v  *  eqn.  (3.  3  5), 

o  3  ' 

and  substituting  this  relation  into  eqn.  (3.  30d).  The  resulting^ equation  may¬ 
be  integrated  once  without  too  much  difficulty,  to  give. 


.41) 


49 


u+  +  2(tanhr+)u*  =  ~~  ((u  V  -  3(u^+)2u  +) 
3r  J  U10  C  1  L 
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(v  +  1) 


8mi4  .r  + 

— ^ —  J  t  tanh  t  dt 
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m  i  3/2  M10  o  J 


2r  lr 


(1^3)  ,„M 

10 


a  «*'  + 

-  1 — —  In  cosh  r 

u10  (^+1> 


(V  +  1)  ,  +  >2 

iJ -  U  (u  )  +  2 

4  10  ,  +;  g8 

lr 


(3.42) 


where  m  ^  and  ju  are  found  from  eqn.  (3.20),  with  the  terms  in  the  expansion 

for  T  calculated  using  eqn.  (3.26).  Another  integration  of  eqn.  (3.42),  to 

find  u*,  would  be  very  difficult  if  not  impossible.  Fortunately,  the  desired 
J  t 

information  can  be  obtained  by  using  eqn.  (3.42)  .  First,  the  expressions  to 

which  u*  must  match  as  r+  +  oo  are  found  from  eqn.  (3.2  3a).  Thus, 

u+(r  +  -  oo,  y)  and  u*(r  +  —  +  oo,  y)  are  known.  If  they  are  used  to  evaluate 

4-  + 

the  left  hand  side  of  eqn.  (3.42),  and  the  known  inner  solutions  (u^  and  ) 
are  used  to  evaluate  the  right  hand  side,  in  the  limits  as  r+  £  oo,  the  re¬ 
sulting  equations  can  be  used  to  obtain  a  relation  for  x^  That  is,  from 


^  The  author  is  grateful  to  Professor  Messiter  for  having  suggested  the  pro¬ 
cedure  used. 
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eqn.  (3,2  3a)  it  is  seen  that 


(ut)  ,  =  A  rxV  +  B  x+  +  C 
3  u,  d  u,  d  u,  d  u,  d 


where  the  subscripts  u  and  d  refer  to  the  limits  r  = 
respectively,  and  where 


=  -  oo  and  r  =  +  oo 


In  addition, 


a  m 
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u  2  lxxo 
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u  il  lxxo  2x  u 
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XXX 
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XXX 
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+  n  *  (0,  y) 


(3.43) 


(3.44a) 


(3.44b) 


(3.  44c) 


(3.  44d) 


(3.45a) 


(3.45b) 


(3.  45c) 
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In  eqns.  (3.44)  and  (3.45),  the  subscripts  u  and  d  used  with  outer  solutions 
refer  to  values  immediately  upstream  and  downstream  of  the  shock  wave 
respectively. 

If  eqn.  (3.42)  is  evaluated  at  r+  -►  -  oo  and  r+  oo,  as  described  above, 
and  coefficients  of  like  powers  of  x+  (i.e.  r+)  on  either  side  of  the  equation 
are  equated,  one  obtains  six  equations,  three  for  each  limit. 

The  coefficients  of  (x+)^  give  equations  which  are  identities,  and  the 
coefficients  of  x+  give  equations  which  may  be  used  to  evaluate  g'  ^(y).  It 
is  the  constant  terms  which  involve  x^  they  are  found  to  be: 
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8  (y  +  1)  u 
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where  the  subscript  i  =  u  when  the  upper  signs  are  used  and  i  =  d  when  the 
lower  signs  are  used.  When  one  equation  is  subtracted  from  the  other,  then, 
one  obtains 


B  ,)  -  (C  + 
d  u 


Cd) 


mf' 


-  1 


(y  +  l)u 


10 


8a  * 


2(*lt  3  2 \ 

(y  +  1) 


(3.47) 


The  desired  relation  is  found  by  substituting  for  B^,  B^,  C^,  and  from 
eqns.  (3.44)  and  (3.45),  using  the  known  outer  solutions  (eqns.  (3.19))  to 
evaluate  (u^  )^,  u^,  etc#  ^ter  considerable  algebraic  manipulation,  it  is 
found  that 
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2(X'<3/2) 

(\  +D 


=  0  (3,48) 


where  is  given  in  eqn.  (3.19e),  in  (3.  37a),  and  ?  in  eqns,  (3.  3e)  and 
(3.  3f). 

At  this  point,  in  the  second  order  calculation,  the  integral  condition 

* 

on  ?  (0,  y)  was  used  to  derive  an  equation  from  which  x  could  be  found. 

*  10 

In  this  third  order  calculation,  the  same  procedure  is  followed,  the  integral 
conditions  being  represented  by  1^,  where 


I ,  =  J  r)**( 0,  y)dy  (3.49) 

0  x 

* 

In  general,  the  governing  equation  for  rj  will  be  a  Poisson  equation,  with 
any  forcing  function  depending  upon  lower  order,  known  functions.  Hence, 
although  the  solution  for  1^  is  not  given  here,  it  is  clear  that  it  can  be  cal¬ 
culated  without  the  necessity  of  performing  any  more  calculations  in  the 
outer  or  inner  shock  structure  regions.  If  eqn,  (3.48)  is  integrated,  then, 
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is  independent  of  Reynolds  number  through  terms  0(e),  and 
thus  is  the  same  as  the  location  of  the  shock  wave  as  Re  oo, 
to  this  order. 

(ii)  Although  the  Hugoniot  jump  conditions  across  a  shock  wave  do 

not  hold  to  the  order  needed,  for  Re  =  0(e  and  therefore 

k  =  O(e),  it  is  possible  to  derive  a  solution  for  x  for  this 
S  £  1 

thick  shock  case. 

(iii)  Although  it  appears  that  x^  depends  upon  m  and  thus  upon 
Reynolds  number  for  k^  =  0(e),  this  is  not  necessarily  the 
case.  That  is,  in  all  probability,  1^  depends  upon  m  also; 
this  dependence  could  be  such  that  x^  is  independent  of  m. 
Support  for  this  possibility  is  given  by  the  fact  that  to  lowest 
order,  dx^  /dy  is  independent  of  m  and  thus  dx^  /dy  =  dx^/dy. 
to  this  order.  If,  indeed,  x^  is  independent  of  m  through  terms 
O(e),  the  location  of  the  sonic  line  in  the  ’’thick11  shock  waves 
found  in  numerical  solutions  in  channel  flows  could  be  used  to 
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locate  the  shock  wave  in  most  practical  cases  even  with  the 
inclusion  of  artificial  viscosity.  In  addition,  if  this  result 
is  indeed  found,  extension  to  more  general  flows  and  high 
flow  Mach  numbers  should  be  considered. 

Finally,  it  should  be  noted  that  composite  solutions  can  be  formed  in¬ 
cluding  now  the  solutions  valid  in  the  inner  structure  region.  In  doing  so, 
it  would  be  necessary  to  write  the  solutions  in  two  parts,  one  for  x  <  x^ 
and  one  for  x  >  x^  . 

Shock  Wave  in  the  Neighborhood  of  the  Throat 

In  the  work  described  in  the  previous  sections,  the  shock  wave  was  lo¬ 
cated  in  the  channel  a  distance  0(1)  downstream  of  the  nozzle  throat.  The 
question  arises  as  to  the  validity  of  these  solutions  as  a  sonic  nozzle  throat 
is  approached  and  the  flow  Mach  number  tends  to  unity.  Because  the  wall 
shape  at  the  throat  is  such  that,  in  general,  f (0 )  =  V  (0)  =  0,  then  in  the 

neighborhood  of  the  throat 

2 

f(x)  =  f(o)  +  ... 

i.  e.,  the  wall  shape  is  parabolic.  Hence,  in  investigating  the  throat  region, 
there  is  no  loss  of  generality  in  considering  a  parabolic  wall  such  that  the 
channel  flow  solutions  are, 

f(x)  =  bx2  (3.51a) 


f(x)  =  +  Bx; 


(3.  51b,  c) 


v2  =  2bxy  (3.  5  Id) 

where  u^,  with  =  0  because  the  flow  is  sonic  at  the  throat,  is  given  by 
eqn.  (3.3a)  and  v^  by  eqn.  (3.  19b).  With  the  above  relations  for  and  v^ 

u  =  1+  eBx  +  ...  v  =  e22bxy  +  ...  (3.  52a,  b) 

and  so  if  x  =  0(a),  say,  in  the  neighborhood  of  the  throat,  then,  expansions 
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for  u  and  v  would  be  such  that  u  =  1  +  O(eof)  and  v  =  0(  €  or).  Because, 
from  equation  (16),  it  is  seen  that 

a2  =  1  -  (u-l)+...  (3.53) 

so  that  U2  -  a^  =  (v  +  l)(u-l)/2,  the  important  terms  in  the  gasdynamic  equa¬ 
tion  (3.  13a)  are  found  to  have  orders  as  shown  in  the  following,  where  the 
order  of  each  term  is  shown  beneath  it: 


,  2  2.  2 
(u  -  a  )ux  -  a  vy 


[vlR(u)  -  (\-l)  R(H)  ] 


€a  •  €a 
a 


2 

€  a 


1 

Re 


2 

a 


Thus,  the  terms  on  the  left  hand  side  of  the  equation,  those  which  are  asso¬ 
ciated  with  channel  flow  solutions,  are  Ofe^c*).  The  viscous  term  on  the 
right  hand  side,  generally  important  in  lowest  order  solutions  only  within 

p*  - 1  /2 

a  shock  structure  region,  is  0(€/Rea).  As  long  as  a  »  (Re  e)  ,  the 

viscous  term  is  negligible,  and  the  channel  flow  solutions  hold  in  lower 

order;  the  viscous  term  could  only  become  important  in  a  thin  shock  structure 

region  imbedded  in  the  region  in  question,  so  the  solution  essentially  remains 

~  -i/2 

that  illustrated  in  the  previous  sections.  However,  when  a  =  0[(Re  e)  ]  , 
then  all  three  terms  are  of  the  same  order.  The  solution  is  not  that  given 
previously,  and  it  is  this  case  we  wish  to  investigate. 

The  region  under  study  in  the  neighborhood  of  the  throat  is  thus  of 
order  a  =  0(Re  e)  7  ]  in  thickness  and  so  we  define  an  inner  variable  and 
velocity  expansions  as  follows: 


ft  =  - 

Or 


u  =  1  +  eat  u  .  (ft ,  y)  +  . . . 


v  =  e 


'a  Vj(ft  ,  y)  +  .  .  . 


(3.  54a) 

(3.  54b,  cl 


2  2  ~  -l 

Because  M  -  1  =  0(u  -  1),  and  a  can  be  written  as  4  =  0[(Re  eu)  J,  it 

is  clear  that 
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«  *  of-  1 1 

[Re(M  -1). 


and  thus  that  a  is  a  measure  of  the  thickness  of  the  shock,  which  fills  the 

whole  region  here,  just  as  k  measures  the  shock  thickness  in  the  previous 

s 

analysis.  Hence,  we  define  a  in  the  same  general  form  as  k  ; 

s 


1 


Q,  =  ~ 

Re  ear  (\  +1) 


1  + 


t) 


(3.  55) 


If  eqns.  (3.  54)  and  (3.  55)  are  substituted  in  eqns.  (3.  10)  and  then  in  eqns. 
(3.9c),  (3.  13c),  (3.  13b)  and  (3.  13a),  it  is  found  that 


R(u)  =  0(e/ar)  =  R(H) 
H  = 


R(v)  =  O(e) 


(3.  56a,  b,  c) 
(3.  56 d) 


-v  2  2 

As  =  0(e  or  ) 


n 


0(e2a2) 


(3.  56  e) 
(3.  56 f) 


and  finally  that  the  governing  equation  for  is, 


M 


1  o,  =  a. 


(3.  57) 


iui£  *  (V+  l)  viy  “  uii& 

Thus,  the  lowest  order  governing  equation  is  the  viscous  transonic  equa - 

U0) 

tion. 

If  S2  "  v  -  u  is  evaluated,  using  eqns*  (3.  56 f)  and  (3*  54),  the 
x  y 

result  is,  to  lowest  order  in  each  term, 

2^1  8“l  2  2 

e  -  +  ...  -  eQ  - —  +  ...  =  0(e  a  )  (3.58) 


^^Cole,  J.  ,  Problems  in  Transonic  Flow,  Ph.  D.  Thesis,  California 
Institute  of  Technology,  1949. 
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Thus,  from  eqn.  (3.55),  the  orders  of  the  two  terms  are  and  Je/ Re 

3/2  ~  -l/2 

If  the  Reynolds  number  is  large  enough,  that  €  »  Re  ,  then 

^lx  =  0°r^i  =  ^  ^(y),  and  such  a  solution  cannot  satisfy  the  matching 

conditions  given  by  the  outer  solutions  for  the  v  velocity  component,  eqn. 

~  - 1/2  2 

(3.5 Id) .  Moreover,  for  this  case,  a  =  (Ree)  «  €  .  Now,  if  the 
outer  solution  for  u  is  expanded  for  x  «  1  to  second  order,  again  for 

wall  shape  and  u^  as  given  in  eqns.  (3.  51)  and  with  u^  as  given  in  eqns. 
(3.  19a)  and  (3.  19c),  we  obtain 


u  ~  1  +  €  bx  + 


+  .  •  . 


Thus,  as  x  0,  the  second  order  term  becomes  the  same  order  as  the 
first  order  term  for  x  =  0(e  *^)  and  x  =  O(c),  indicating  the  possi¬ 

bility  that  solutions  in  two  more  regions  should  be  considered.  However, 
it  was  shown  in  reference  2  that  it  Ls  possible  to  write  the  solution  for  u 
in  a  form  such  that  the  singular  terms  in  u^  do  not  arise,  being  contained 

in  u.  through  the  use  of  an  expansion  for  C  .  Hence,  solutions  in  the 

1  1/2  W 

region  0(6  )  in  thickness  do  not  differ  significantly  from  the  channel 

flow  solutions  already  being  used.  However,  in  the  thinner  region 
(  x  =  O(e)),  it  is  clear  that  the  flow  is  two  dimensional  even  in  lowest 

order  and  thus  quite  different  from  the  outer  channel  flow  described  by 
eqns.  (3.  19a)  and  (3.  19c);  the  governing  equation  for  inviscid  flow  is  the 
nonlinear  transonic  small  disturbance  equation.  Hence,  for  a  «  €, 

the  throat  region  under  consideration  would  be  contained  within  the  region 
of  order  c  in  thickness  and  the  M  outer1'  solutions  to  be  used  for  match¬ 
ing  would  be  the  solutions  in  the  latter  region,  not  the  outer  channel 
solutions  given  in  the  previous  section. 

If  a  =  O(c),  then  v^a  =  and  a  velocity  potential  may  be  intro¬ 

duced  into  the  viscous  transonic  equation  (3.  57);  the  lowest  order  solu¬ 
tions  for  u  are  two  dimensional  and  the  outer  solutions  to  which  the  inner 
solutions  must  match  are  those  given  in  eqns.  (3.52a,  b). 
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0,  so 


Finally,  if  a  »  e,  then  from  equation  (3.58),  u  ^  = 

=  &j(&).  In  this  case,  for  which  Re  «  e  ^  is  the  condition,  the  ve¬ 
locity  components  must  match  with  the  outer  solutions,  eqns.  (3.  52a) 

and  (3.52b)  as  |x|  0  and  |x  |  -►  oo.  The  Reynolds  number  under  con- 

~  -2  1/2 

sideration.  Re  =  0(e  )  so  that  a  =  0(e  ),  fulfills  this  condition  and 

so  it  is  this  case  which  will  be  analyzed  here.  It  will  be  shown  that  it 
is  possible  to  derive  an  exact  solution  to  the  viscous  transonic  equation 
for  an  arbitrary  (parabolic)  wall  shape. 

It  may  be  noted  that  because  the  whole  idea  of  these  calculations 
is  to  derive  solutions  comparable  to  the  thick  shock  wave  results  found 
in  numerical  solutions,  it  is  necessary  to  consider  a  Reynolds  number 
small  compared  to  those  associated  with  typical  channel  flows.  Never¬ 
theless,  the  solutions  to  be  shown  are  of  more  general  interest  than  the 
limited  Reynolds  number  range  would  indicate.  That  is,  they  may  be 
used  to  illustrate  the  manner  in  which  a  shock  wave  forms  at  the  throat 


and  then  evolves  into  a  thin  wave  as  it  moves  downstream  in  response  to 
a  pressure  condition  imposed  on  the  flow  downstream  of  the  shock  wave. 
Thus,  the  solutions  will  match  not  only  with  the  outer  channel  flow  solu¬ 
tions  when  the  shock  wave  is  within  this  inner  throat  region,  but  also 
with  the  solutions  including  the  shock  structure  when  the  shock  wave  has 
moved  downstream  of  the  throat  into  the  overlap  region. 

For  the  case  considered,  then.  Re  =  0(e  and  from  eqn.  (3.55), 

~  -2 

if  the  same  constant  of  proportionality  between  Re  and  e  is  to  be  used 
as  in  the  previous  section  (eqns.  ( 3.  8)  with  kg  =  me), 

a  =  (me)1/2  (3.59) 

Now,  from  eqn.  (3.58),  u  ^  =  dj(x)  and  so  from  eqn.  (3.  57^, 

*iy  =  30 

0  1  =  Fj(ft)y  +  F2$)  (3.60) 

Also,  as  |x  |  —  oo,  v  =  €2<*v^  +  ...  (with  as  in  eqn.  (3.60))  must 
match  with  eqn.  (3.  52b),  written  in  terms  of  inner  avariables.  Hence, 

F2(fc)  =  0  and  F^ft)  =  2b£  so  that  apparently 
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v  j  =  2bx  y  (3,6  1) 

throughout  this  inner  region,  subject  to  the  boundary  conditions  at  the  wall 

being  satisfied.  These  conditions  are  given  by  eqn.  (3.  17)  with 
2  2  2 

y  =  1  +  e  f(x)  =  1  t  e  bx  ,  where  eqn.  (3.  51a)  has  been  used  for  f(x). 

w 

Thus,  the  boundary  conditions  at  the  wall  reduce  to 

0  (ft,  +  1)  =  +  2b&  (3.62) 

and  are  indeed  satisfied  by  eqn.  (3.6  1). 

As  noted  previously,  y^  is  taken  to  include  the  displacement  thickness 
of  the  boundary  layer.  Again,  the  effects  of  the  interaction  between  the  shock 
wave  and  the  boundary  layer  may  be  estimated  using  the  solutions  valid  within 

(9) 

the  boundary  layer  in  an  interaction  region  and  the  orders  of  the  velocity 

perturbations  (ea  for  the  u  component  and  e^a  for  the  v  component)  in  the 

flow  outside  the  boundary  layer;  again,  it  is  easily  shown  that  these  effects 

are  small  compared  to  the  change  in  displacement  thickness  of  the  undisturbed 

2 

boundary  layer  over  the  interaction  region.  Since  e  a  v  the  lowest  order 
solution  for  v  within  this  thick  shock  region,  satisfies  the  boundary  condition 
associated  with  the  undisturbed  boundary  layer  displacement  thickness,  it 
is  only  in  higher  order  terms  that  corrections  might  be  needed  (e.  g.  ,  terms 
O(e^cr^)  =  O(e^)),  because  these  boundary  conditions  were  not  satisfied. 

Such  higher  order  terms  are  not  considered  here. 

With  v  j  as  given  in  eqn.  (3.6  1),  the  governing  equations  for  u  eqn. 

(3.  57),  can  be  written  as  follows: 


a  A  2  bx  A 

u  .  u  a  -  : — 777  =  u  .  aa 
1  lx  (\  +1)  1 XX 


which  can  be  integrated  once  to  give 


i$ 


A?- 

U1 


u2  a  2 
B  x 


2BA 


(3.63) 


where  the  constant  of  integration  has  been  written  as  -2 BA  for  later 
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convenience.  Equation  (3.6  3)  can  be  transformed^^  into  a  standard  form  of 
Weber' s  equation,  the  solution  of  which  can  be  written  in  terms  of  parabolic 
cylinder  functions^*2).  In  terms  of  the  standard  solutions,  U(A,  X)  and  V(A,  X) 
given  in  reference  12,  where 

X  =  \/~B  x,  (3.64) 

the  solution  for  Uj  maybe  written  in  two  parts.  For  x  >  0, 


ux  X[U (A,  X)  -  CjVf^X)]  +  (2A+1)U(A+1,X)  -  C  ^2 A-l) V(A- 1,  X) 

7b  =  U(A,  X)  +  C1V(A,X) 

and  for  x  <  0 


1  .  v  „A.n  U(A+1.  -X) 

Tb  X  '  *2A+1)  U  (A,  -X) 


(3.66) 


where  in  order  that  the  solutions  agree  at  X  =  0, 


ci  = 


Zr(j  -  A) tan  [*(^+^)] 

1  -  tan2[>  (  “  +  ~  )] 

(12) 


(3.67) 


In  eqn.  (3.67),  T(z)  is  the  Gamma  Function 

Equation  (3.6  3)  was  obtained  and  solved  first  by  Kopystynski  and 
(13) 

Szaniawski  in  their  study  of  flow  in  a  nozzle  throat.  They  wrote  the 
solution  in  terms  of  confluent  hypergeometric  functions;  the  solution  here, 
written  in  terms  of  parabolic  cylinder  functions  is  equivalent  to  theirs  and 
somewhat  more  convenient  to  use.  It  is  interesting  to  note  that  Kopystynski 
and  Szaniawski  showed  velocity  distributions  which  appeared  to  be  those  for 


Murphy,  G.  M. ,  Ordinary  Differential  Equations  and  Their  Solutions, 
Van  Nostrand  Company,  Inc.,  1960. 

^^Handbook  of  Mathematical  Functions,  Eds.  M.  Abramowitz  and  I.  A. 
Stegun,  N.B.S.  Applied  Mathematics  Series  .55,  1964. 

^^Kopystynski,  J.  and  Szaniawski,  A.,  Structure  of  Flow  in  a  Nozzle 
Throat,  Archiwum  Mechaniki  Stosowanej,  3,  17  (1965),  453-466. 
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very  thick  shock  waves  but  did  not  evaluate  their  solution  in  the  limit  as 
this  relatively  rapid  variation  in  velocity  moves  downstream  of  the  throat. 
Thus,  they  did  not  show  that  this  solution  evolves  into  that  for  a  thin  shock 
wave  across  which  Hugoniot  jump  conditions  hold.  It  is  this  point  which 
is  the  subject  of  the  present  section. 

It  is  possible,  then  to  write  an  exact  solution  for  u,  in  terms  of  known 
functions.  Now,  eqn.  (3.6  3)  is  a  relatively  simple  first  order  nonlinear  dif¬ 
ferential  equation  which  it  is  quite  easy  to  integrate  numerically  using  only  a 
programmable  hand  calculator.  Nevertheless,  the  analytical  solution  is  use¬ 
ful  in  that  it  may  be  used  to  demonstrate  that  it  becomes  exactly  that  given 
in  the  previous  section  as  the  shock  wave  moves  downstream  of  the  throat. 

The  expansions  necessary  for  matching  with  the  outer,  channel  flow 
solutions,  are  found  by  writing  U  and  V  for  |X|  »  1.  From  reference  12, 

one  finds  that  for  X  large  and  A  moderate,  for  X  »  |Aj, 


U  (A,  X)-  e 


-A 


X 


1 

2 


(A+  (A+  ~) 


2X 


V(A,  X) 


(A-  ~)  (A-  I*) 


2X 


Hence,  from  eqns.  (3.64),  (3.65),  and  (3.66),  as  x  -►  oo 


-  Bx 


(2  A- 1) 

T  .  .  • 

A 

x 


(3 


(3 


(3 


and  as 


$ 


00 


&  =  b£  +  +  . . . 

1  A 

X 

Equations  (3.69)  and  (3.70)  can  be  checked  by  substitution  into  eqn.  (3.6  3) 
and,  indeed,  could  have  been  derived  from  it.  Now,  the  present  solutions 
hold  for  the  case  where  the  shock  wave  is  in  the  neighborhood  of  the  throat. 


.68a) 


.68b) 


.69) 


.70) 
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Although  the  flow  upstream  of  the  throat  could,  in  general,  be  either  subsonic 
or  supersonic,  we  consider  here  the  case  where  it  is  subsonic.  Hence,  as 
&  -  oo  and  &  -*■  +  oo,  the  inner  solutions  should  match  with  outer  channel 

flow  solutions  which  are  valid  for  subsonic  flow  upstream  of  the  throat  and 
subsonic  flow  downstream  of  the  shock  wave,  respectively,  both  written  in 
the  limit  as  x  0.  For  the  solutions  valid  downstream  of  the  shock  wave, 
this  means  also  that  x  0  and  as  a  result,  u,^  -►  0.  The  required  outer 
solutions  are  found  from  eqns.  (3.  18a),  (3.  19a),  (3.  54a)  and  (3.  59),  using 
eqns.  (3.  51)  for  f(x),  and  u^.  Thus,  for  x  -*■  0  from  upstream  of  the  throat 


C  +  mB 

u  =  1  +  ea  <  B&  +  - 

mBx 


and  for  x  0  from  downstream  of  the  throat. 


mBx 


+  .  •  • 


(3-71) 


u  =  1  -  €<*  \  Bx  + 


(3.72) 


where  both  equations  are  written  in  terms  of  the  inner  variable  x  .  1/latching 
the  inner  and  outer  solutions,  one  finds  that, 


C  =  2mBA  =  C,,  (3.73) 

2u  2d 

If  this  result  is  compared  with  that  given  in  eqn.  (3.  39),  it  is  seen  that  for 

u.  _  =  Bx  =  a  Bx  =  (me)^^Bx  ,  the  two  results  agree,  in  lowest 

10  so  so  so 

order.  Thus,  in  the  throat  region,  choosing  a  value  for  C  =  C  is 

2d  2u 

equivalent  to  choosing  a  value  for  A,  and  as  will  be  seen,  this  means  setting 
the  location  of  the  shock  wave  just  as  in  the  outer  region. 

For  A  =  -  1/2,  Cj  =  0  (eqn.  3.67)  and  the  solutions  for  u  ^  (eqns. 

(3.6  5)  and  (3.66))  reduce  to 


u  =  n/B  X  =  Bx 


(3.74) 


This  is  the  solution  for  a  flow  accelerating  from  subsonic  to  supersonic  ve¬ 
locities  with  no  shock  wave;  that  is,  any  shock  wave  which  might  occur  is 
positioned  downstream  of  the  throat  region  under  consideration.  On  the  other 
hand,  for  values  of  A  greater  than  -  1/2,  solutions  for  u  ^  show  shocklike 


behavior,  as  illustrated  in  figure  2.  The  solutions  in  figure  2  were  found  by 
integrating  eqn.  (3.6  3)  numerically,  with  an  H.  P,  29C  programmable  cal¬ 
culator,  using  eqn.  (3.70)  to  obtain  initial  values  for  u  ^/VB.  They  could 
have  been  obtained  also  by  using  eqn.  (3.6  5)  and  (3.66)  with  numerical  values 
for  U(A,  X)  and  V(a>X)  obtained  from  tables  in  reference  12.  The  solutions 
show  that  there  is  a  certain  value  for  A,  (somewhat  larger  than  -  0.  30) 
greater  than  which  there  is  no  shock  wave  since  no  supersonic  flow  exists. 
The  flow  is  simply  a  viscous  channel  flow.  As  A  decreases  below  this  value 
and  tends  toward  -  1/2,  a  shock  wave  forms  and  moves  downstream.  The 
more  closely  A  approaches  -  1/2,  the  closer  is  the  flow  upstream  of  the 
jump  in  velocity  to  the  supersonic  inviscid  channel  flow  solution 
u^  =  =  B&  and  the  thinner  is  the  region  associated  with  the  jump  in 

velocity.  Downstream  of  the  jump,  the  solution  approaches  the  subsonic 
inviscid  channel  flow  solution,  u  ^  =  -  \fB  X  =  -  Bx  .  The  fact  that  the 

Hugoniot  shock  wave  jump  conditions  are  not  satisfied  in  the  immediate 
region  of  the  jump  is  clear;  i.  e. ,  these  jump  conditions  are  u  ^  =  -  u 
and  in  figure  2  it  can  be  seen  that  the  solutions  are  not  symmetric  about 


The  solutions  shown  in  figure  2  apparently  show  the  evolution  of  a  shock 
wave  as  it  forms  at  the  throat  and  moves  downstream  in  response  to  a  pres¬ 
sure  condition  impressed  upon  the  flow  downstream  of  the  shock  wave. 

In  order  to  prove  this  contention,  it  must  be  shown  that  as  the  shock  wave 
moves  downstream,  this  inner  solution  matches^  with  the  solutions  valid 
for  the  case  where  the  shock  wave  is  located  at  a  distance  of  0(1)  downstream 
of  the  throat,  but  written  now  in  the  limit  as  the  shock  wave  moves  toward 
the  throat.  From  the  solutions  shown  in  figure  2,  it  is  clear  that  the  inner 
limit  process  involves  A  -*■  -  1/2  as  €  0. 

The  outer  solutions  to  which  the  inner  solution  must  match  is  a  com¬ 
posite  solution,  valid  to  first  order,  formed  from  the  solutions  found  in  the 
previous  section.  Thus,  if  the  flows  upstream  of  the  shock  wave  and  in  the 
inner  structure  regions  are  considered,  the  composite  solution  is  the  sum 
of  the  channel  flow  and  inner  structure  solutions  minus  the  common  (found 
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from  matching)  terms. 


u  =  1  +  e 


u^x)  +  ei^  (x  )  -  €Uj0  +  ... 


(3.75) 


Equation  (3.75)  is  valid  for  x  <  x^  ,  (x  <  0),  because  the  supersonic  value 

of  u.(x)  was  used  in  determining  the  common  term,  '"’he  corresponding  equa- 

*  +  + 

tion  valid  for  x  >  x  would  be  u  =  1  +  eu,(x)  +  eu  (x  )  +  eu  +  .  .  .  .  Now, 

-I-  *  h  ^ 

because  x  =  (x-x.  )/k  and  x*  «  x,  +  ex.  +  €  *,,+•••#  with  k  =  O(e), 

is  i  io  il  i 1  s 

it  is  seen  that  for  eqn.  (3.7  5)  to  be  uniformly  valid  to  O(e),  it  is  necessary 

only  to  include  the  first  two  terms  in  the  expansion  for  x^  ,  i.  e. , 

x,  =  x#  +  ex.  _,  and  this  is  understood  to  be  the  case  in  what  follows. 

Finally,  because  when  the  shock  wave  is  within  the  throat  region,  matching 
between  the  inner  throat  region  solution  and  the  outer  channel  flow  solutions 
has  already  been  demonstrated,  it  only  remains  to  demonstrate  that  the 
solutions  match  in  the  immediate  vicinity  of  the  shock  wave  as  it  moves 
toward  the  throat.  Therefore,  the  outer  solution,  eqn.  (3.75)  is  written  in 
the  limit  as  x  -►  x^  ,  for  x^  «  1  such  that  x**  =  0(1).  Then,  to  lowest  order, 
eqn.  (3.7  5),  and  uin  become 


10  + 

u  =  1  -  €  u^0  tanh  (~j—  x  )  +  .  . 


(3.76a) 


(3.  76b) 


where  for  x  «  1,  =  Bx,  and  where  the  solutions  for  x^  ,  eqn.  (3.  33a)  has 

been  employed. 

The  corresponding  limit  to  be  used  in  evaluating  the  inner  solution  is 

X  =  Xi  +X  -  Xi  (3.77a) 


X  -  X. 


«  1  «  X , 


A  2  2 


u  =  w(e)  «  1 


(3.77b) 


(3.77c) 


( 3.  77d) 
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When  eqn.  (3.67),  is  evaluated  using  eqn.  (3.77c),  for  A,  it  is  found  that 


C 


1 


W  If 


+  0(w2) 


(3.78) 


Next,  if  eqn.  (3.78)  and  the  asymptotic  expressions  for  U(A,  X)  and  V(A,  X), 
eqns.  (3.68a)  and  (3.68b)  respectively,  are  substituted  into  the  equation  for 
u^  valid  for  x  >  0,  eqn.  (3.65),  the  following  result  is  obtained  for  X  »  1, 
and  thus  x  »  1. 


-  X  * 


u'flv)  _1X  e  2 


1  -  (oN2t) 


-  X  + 


l+lu^l’^e  2  +X"2+... 


(3.79) 


Now,  starting  with  eqn.  (3.77a)  and  using  the  definitions  of  X,  x  and  x+, 
eqns.  (3.64),  (3.54a)  and  (3.21a)  respectively,  one  can  write  X  as 


X 


m 


x^  ^  +  *J~€ mB 


(3.  80) 


where  x^  has  been  expanded  using  only  the  first  two  terms,  as  mentioned 
previously.  In  addition,  eqn.  (3.76b)  has  been  used  for  and  eqn.  (3.  59) 
for  a  .  When  eqn.  (3.  80)  is  substituted  into  eqn.  (3.  79),  it  is  seen  that  if 
to  lowest  order 


then. 


,  .  10 

“(e)  =  77rnI7c  “P 

again  to  lowest  order. 


(3.81) 


+ 

-  u10  tanh  r  (3-  82) 

where,  again,  rT  =  u^xx/2.  Therefore,  this  solution  from  the  inner  throat 
region  matches  the  corresponding  outer  solution  in  the  immediate  vicinity  of 
the  shock  wave,  eqn.  (3,76a).  Evidently,  the  solutions  found  in  this  section 
do  represent  the  evolution  of  a  shock  wave  as  it  moves  downstream  of  the  throat. 
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changing  from  a  thick  shock  across  which  Hugoniot  conditions  are  not  satisfied 
to  a  thin  shock  wave  across  which  they  are  satisfied. 

Finally,  it  is  interesting  to  note  that  the  solutions  shown  in  figure  2  are 

(14) 

quite  similar  to  those  given  previously  for  steady  channel  flow  by  Sichel 
(two-dimensional)  and  Sichel  and  Yin^^  (axis ymme trie),  and  for  unsteady 
two-dimensional  flow  by  Adamson  and  Richey  .  In  all  of  these  analyses, 
a  similarity  transformation  was  used  to  reduce  the  viscous  transonic  equation 
to  an  ordinary  differential  equation.  Then,  in  references  14  and  15,  numeri¬ 
cal  solutions  were  obtained;  in  reference  16,  numerical  results  were  obtained 
for  thick  shock  waves  and  analytical  solutions  for  thin  shockwaves.  In  all 
cases,  as  the  solution  for  the  perturbation  in  u  went  through  the  jump  caused 
by  the  shock  wave,  it  overshot  the  solution  associated  with  decelerating  sub¬ 
sonic  flow  downstream  of  the  shock  and  then  approached  this  solution  asymp¬ 
totically  from  beneath.  As  seen  in  figure  2,  the  present  solutions  do  not  show 
this  feature;  instead,  u^  varies  monotonically  from  its  peak  value  to  the 
solution  for  subsonic  decelerating  flow  (u^  -  -  Bx  ).  This  difference  occurs 
because  the  similarity  solutions  for  simple  supersonic  accelerating  flow  and 
subsonic  decelerating  flow  are  not  symmetric  about  u^  =  0.  That  is,  at  any 
given  location  the  solution  for  subsonic  decelerating  flow  is  not  the  negative 
of  the  solution  for  accelerating  supersonic  flow.  Hence,  these  solutions  can¬ 
not,  in  themselves,  satisfy  the  jump  condition  across  a  shock  wave, 

u,  ,  =  -  Q,  ;  an  overshoot  in  the  solution  results  when  a  shock  wave  occurs. 

Id  lu 

(14) 

Sichel,  M. ,  The  Effect  of  Longitudinal  Viscosity  on  the  Flow  at  a  Nozzle 
Throat,  J.  Fluid  Mech. ,  25  (1966),  769-786. 

^^Sichel,  M.  and  Yin,  Y.  K. ,  An  Axisymmetric  Similarity  Solution  for  Vis¬ 
cous  Transonic  Nozzle  Flow,  J.  Fluid  Mech.,  28^  (19 67) ,  512-522. 

Adamson,  T.  C.,  Jr.  and  Richey,  G.  K. ,  Unsteady  Transonic  Flows  with 
Shock  Waves  in  Two-Dimensional  Channels,  J.  Fluid  Mech.,  60_  (1973), 
363-382. 
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Perhaps  the  most  important  difference  between  the  similarity  solutions 
and  those  presented  here  is  that  the  former  cannot  be  employed  for  arbitrary 
wall  shapes.  Thus,  a  wall  must  be  associated  with  a  streamline  calculated 
from  the  similarity  solution;  no  arbitrariness  is  allowed.  Not  only  does 
this  limit  the  usefulness  of  similarity  solutions  in  practical  applications,  it 
leads  to  wall  shapes  with  bends  or  angles  at  the  point  where  the  streamline 
associated  with  the  wall  passes  through  the  shock  wave,  depending  upon  the 
thickness  of  the  shockwave.  The  present  solutions  allow  arbitrary  wall 
shapes  to  be  considered;  in  those  cases  where  the  streamline  variation 
through  the  shock  wave  is  different  from  that  of  the  wall  (including  the  ef¬ 
fects  of  the  displacement  of  the  boundary  layer),  i.  e.,  when  interaction  ef¬ 
fects  are  important,  corrections  can  be  added. 
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